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FOREWORD 


This report summarizes the work completed on the research project 
“Radiative Transfer Models for the Earth Radiation Budget Studies." The 
work was supported by the NASA/Langley Research Center (Experiment .Analysis 
Branch of the Atmospheric Environmental Sciences Division) through research 
grant NAG-1-21. The grant was monitored by Dr. John T. Suttles of the 
Atmospheric Environmental Sciences Division. The authors are grateful to 
Dr. S. K. Gupta for providing various help during the course of this study. 
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ANALYSIS OF LONGWAVE RADIATION FOR 
THE EARTH-ATMOSPHERE SYSTEM 


By 

S. N. Tiwari, 1 C. S. Vemuru, 2 and S. V. Subramani an 3 


SIMMARY 

Accurate radiative transfer models are used to determine the upwelling 
atmospheric radiance and net radiative flux in the entire longwave spectral 
range. The validity of the quasi-random band model is established by com- 
paring che results of this model with the results of line-by-line furmul a- 
tions and with available theoretical and experimental results. Existing 
radiative transfer models and computer codes are modified to include various 
surface and atmospheric effects (such as surface reflection, nonequilibrium 
radiation, and cloud effects). The program is used to evaluate the radia- 
tive flux in clear atmosphere, provide sensitivity analysis of upwelling 
radiance in presence of clouds, ana determine the effects of various clima- 
tological parameters on the upwelling radiation and anisotropic function. 

Tne program is used a 1 so to evaluate homogeneous and nonhomogeneous gas 
emissivities under different conditions and the results are presented in a 
supplement to this report entitled "Accurate Evaluation of Homogeneous and 
Nonhomogeneous Gas Emissivities." 
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1. INTRODUCTION 


Extensive study of the radiative transfer phenomena in the Earth's at- 
mospheric system has been carried out in the last two decades (refs. 1 to 
4). This is important for the understanding of the meteorological process 
on all scales and the spatial variation in surface temperature in the Earth 
atmosphere. Techniques for measuring the Earth's surface temperature in- 
clude airborne instruments and satellite-mounted radiometers. In order to 
understand and interpret the instrument performance and readings, it is 
desirable to develop radiation models and nunerical techniques that accou.it 
for the absorption and attenuation of actual atmospheric radiation. Devel- 
opment of accurate models for radiative transfer in the atmosphere is ex- 
tremely important for Earth radiation budget studies ana climate modeling 
(refs. 4 to 6). These models have to be used for simulation ana interpreta- 
tion of Earth radiation budget measurements as well as for retrieval of 
various surface and atmospheric parameters from satellite-measured radiances 
(ref. 7). Since the radiation budget of the planet has been identified as 
an important element of the climate system, its measurements are being at- 
tempted with increasing accuracy (refs. 8 and 9). As a result, considerable 
improvement is warranted in the accuracy of the theoretical models dealing 
with atmospheric radiation transfer. 

Many models for radiation absorption by molecular gases are available in 
the literature. The simplest one is the gray gas model (or the emissivity 
approximation) and the most sophisticated and accurate one is the line-by- 
line { LBL ) model (or the direct integration procedure). Between the emis- 
sivity approximation and direct integration method lie several narrow-and 
wi de-band models and band model correlations which vary greatly in complex- 


2 


> 


ity and accuracy. A comprehensive review of various line and band models is 
available in reference 4. Use of either a LBL model or a narrot'-band mod- 
el is suggested for most atmospheric applications. The narrow-band models 
usually recommended for atmospheric studies are the Elsa’sser (or regular 
model, statistical (Meyer-Goody or Goody) model, and quasi-random band ( QRB) 
model. The ORB is probably the best band model to represent accurately the 
absorption of a vibration-rotation band and is suitable for calcula- 
ting the atmospheric transmittance and upwelling radiance. The fundamental 
features of the QRB are discussed, in detail, in references 10 to 12, and 
the procedure for calculating the atmospheric transmittance and upwelling 
radiance is given in reference 12. In spectral ranges where both line ab- 
sorption and scattering are important, a widely used approximation for cal- 
culating spectrally integrated radiative flux is the cxponenti al-sum fitting 
of transmissions (ESFT) method. The basis for this method is that the 
transmission function for a given spectral Interval is fit by a sum of ex- 
ponentials. The method is described in reference 13. 

Radiative transfer models used in earlier climatic investigations em- 
ployed radiation charts, generalilzed absorption coefficients, and emissiv- 
ity approximations (refs. 5, 14-16). Rodgers (ref. 17) has indicated that 
the use of mul ti-interval narrow-band radiative transfer schemes in climate 
modeling studies will constitute a significant step forward and result in 
improved accuracy of the model output. Fels and Kaplan (ref. 18) have in- 
vestigated the effects of using different radiative transfer schemes on the 
thermal structure of the atmosphere and its consequences to atmospheric dy- 
namics. They employed two different radiative absorption models, the emis- 
sivity approximation and Goody's statistical band formulation, and performed 
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numerical experiments with the NCAR general circulation model. They ob- 
served a significant difference in the cooling rates in the two experiments 
which resulted in signif icant ly different mean temperature fields and merid- 
ional circulation. 

Very high accuracy can be achieved in the radiation computation by 
using the LBL integration procedure in the radiative transfer models (ref. 
19). However, the procedure is too cumbersome and makes excessive demands 
on computer time. Tiwari and Gupta (ref. 20) have shown that the QRB model 
can be used for computing atmospheric transmittances with accuracy compara- 
ble to that of the LBL method and with computer usage more than an order of 
magnitude smaller. Kunde (ref. 11) has also used this model to compute out- 
going infrared radiances fro..* planetary atmospheres. However, before use of 
the QRB model can be recommended for Earth radiation budget and climate 
modeling studies, further work needs to be done to val idate the model on a 
sound basis. This model should be used for absorption bands of different 
species in different spectral ranges. It is quite possible that the model 
is not justified at shorter wavelengths and smaller pressure path lengths. 
Furthermore, under realistic atmospheric conditions, the model may give good 
results in certain spectral ranges, but is poor in other range:.. 

For critical applications, it is essential to validate the quasi-random 
band model under as many different but realistic conditions as possible by 
comparing the results of this model with available experimental and theoret- 
ical results. For several molecular species, experimental results for spec- 
tral transmittance and total band absorpance are given by Burch et al. 

(ref. 21) under different pressure and path length conditions. Tn us it is 
highly desirable to compare the results of the QRB model with these exoeri- 
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mental results. For cases where experimental results are not available, it 
is important to compare the QRB results with LBL results. For certain spec- 
tral ranges and atmospheric conditions, results of atmospheric trariamit- 
tances are available in the literature which have been obtained by using a 
sophisticated program called LOWTRAN (ref. 22). It is therefore desirable 
to obtain the QRB results exactly for these conditions for comparison with 
the LOWTRAN results. Another model ised frequ.. it ly in atmospheric studies 
is the K-distribution approximation (refs. 23-27). It has been applied suc- 
cessfully only to water vapor bands (refs. 24 and 25), but there are indica- 
tions thct it could be usad with reasonable accuracy for other bands also 
(refs. 26 and 27). It is therefore important to compare the results of the 
QRB model with the K-distribution formulation for different bands under 
varying conditions. Recently (refs. 28 and 29), another model and computer 
code called FASCODE (Fast Atmospheric Signature Code) has Seen developed f 
the line-by-line calculation of radiance and transmittance with particular 
applicability to the Earth's atmosphere. In this model, an algoritrrn for 
the accelerated convolution of line shape functions (Lorentz, Voigt and 
Doppler) with spectral line data is used. The contribution from continuum 
absorption is also included in the irndel. Tne program is applicable to 
spectral regions from the microwave to the visible. It may, therefore, be 
desirable to compare the results of this model with the results of QRB for- 
mui ation. 

The objective of this study is to validate the QRB model for a few real- 
istic conditions by comparing the results of this model with the LBL, exper- 
imental, and LOWTRAN results. Detailed verification of the QRB model re- 
sults with the results of other formulations is beyond U,e scope of this 
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study. After the model validations, the aim of this study is to use the (JR B 
model for evaluating the gas emissivities in several important spectral 
ranges and also in the entire longwave range. The QRB formulation is espe- 
cially useful for this when a mixture of several molecullar species are in- 
volved. This study could provide benchmark solutions for gas emissivity 
under different pressure and temperature conditions. Such information are 
useful not only in atmospheric studies, but also in the fields of infrared 
signature work, combustion processes, and fire research. Another important 
aim of this study is to use the QRB model for evaluating the upwelling at- 
mospheric radiance under different realistic surface and atmospheric condi- 
tions. Such formulations (and relevant information) are very useful in de- 
veloping an accurate data reduction scheme for the measurement of atmospher- 
ic pollutants by remote sensing (refs. 4 and 30), and for Earth radiation 
budget and climate modeling studies (refs. 5-9, 31-33). 

A major factor influencing the radiation balance and the general circu- 
lation of the Earth's atmosphere is the presence of clouds which occupy 
about 50 percent of i.ne planet Earth on a global basis. Clouds absorb and 
scatter the incoming solar radiation and absorb and emit terrestrial radia- 
tion. Although clouds have been included in the study of transfer of solar 
radiation through the Earth's atmosphere for many years, there have been 
very few studies which include the effects of clouds on the longwave radia- 
tion (refs. 32-40). Clouds vary greatly in thickness, height, liquid water/ 
ice content, and geometrical shape and size; and ail these factors contrib- 
ute in a complicated manner to the large variability of the cloud radiative 
properties. Some of these effects are being considered in recent studies 
(refs. 41-47). Most of the treatments given to clouds in the long-wave 
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radiation transfer models have been very simplistic. It is therefore essen- 
tial to modify the existing radiative transfer models appropriately and in- 
vestigate the effects of large variability in the radiative and geometrical 
properties of the clouds on the thermal radiances and fluxes. 

The basic formulation of the radiative transfer equations and the ex- 
pressions for the upweTMng radiance and flux are presented in the next 
chapter, "Basic Theoretical Formulation." The spectral models used in this 
study are discussed briefly in chapter 3 entitled "Absorption Models." The 
computational procedure and data source for calculating the transmittance, 
upwelling radiance, and radiative flux are given in chapter 4. The results 
of the entire study pertaining to atmospheric applications are presented in 
three separatj chapters entitled, "Accurate Evaluation of Longwave Radiation 
in Clear Atmosphere," "Sensitivity Analysis of Upwe^ing Radiance in Pres- 
ence of Clouds," an^ "Evaluation of Anisotropic Functions in the Longwave 
Region." The results of emissivity calculations are given in a supplement 
to th' c report entitled, "Accurate Evaluation of Homogeneous and Ncnhomogen- 
eous Gas Emissivities. 
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2. 8ASIC THEORETICAL FORMULATIONS 


Basic governing equations for radiative transfer in the atmosphere are 
available in the literature. However, many of these equations need to be 
modified for specific applications. For some applications, entirely new re- 
lations are needed to express a particular phenomena. 

The radiation emergent from the atmosphere may be given by the expres- 
sion (refs. 4, 30): 

E(«) = e g m + E r («) + E gr (») + E, (a) + e 4R («) (2.1) 

The various components of the upwelling radiation are pictori lly shown in 
figure 2.1 and are defined in the list of synbols. In genet al, these are 
functions of surface and atmospheric temperatures, surface emittance and re- 
flectance, sun zenith angle, scattering characteristics of particles, and 
transmittance of the atmosphere. 

In the spectral region of infrared measurements, the effect of scatter- 
ing and solar-reflected radiation is usual ly omitted. Hence, the expression 
for thermal radiation emerging from a plane-parallel atmosphere can be writ- 
ten as 

E(«) = E g M + E gr (<*>) = e(w) B(<u,T s ) t(co,0) 


( 2 . 2 ) 

h 

+ / q B[u,T(z) J[dx(o),z)/dz]dz + p(<i>)F-(o>,T) t(u>,0) 


The first term on the right-hand side of equation (2.2) represents the 
radiation from the surface; the second term is the radiation from the atmos- 
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phere, and the third term represents the reflected component of the downward 
radiation. The contribution of the reflected atmospheric radiation from the 
surface is usually neglected fo* surfaces with relatively high values of 
surface emittance and for the spectral regions where the downward atmos- 
pheric emission is small. 

The contribution of the sunlight reflected from the surface from the 
surface is important at shorter wavelengths and is given by the component 
E r (w) as 

E R (w) * - [l -£(«)] cos 9 H (oj) [t (to) ]“ (2.3) 

ir s 

where 0 is the sun's zenith angle, [l - e(u)j is the ground reflectance of 
the surface, H s (w) is the sun irradiance on the top of the atmosphere, 
a = 1 + f(0) where f(0) = sec 0 for 0 < 0 < 60° and Ch 0 for 0 < 60* with 
Ch 0 denoting Chapman's function, and t(uj) = t(oj,0) is the transmission 
“ertically through the atmosphere. 

For radiation budget and cooling rate calculations, however, the 
required quantity is the flux density. Upward flux density can be obtained 
precisely by integrating the upwelling radiance ever the zenith angle 6 
and the azimuth such that 

2n /2 

F(«,h) = J Q d4> /*'£(<*>, h) sin 0 cos 0 d0 (2.4) 

Integration of equation (2.4) by using detailed angular distribution of 
radiation is a tedious problem. However, it is simplified considerably for 


9 


F 


a plane-parallel atmosphere and assuming that the source function in equa- 
tion (2.2) is isotropic. It is possible with the above assumption to adopt 
the two-stream approximation whereby the equations of transfer are reduced 
to only two. 

For the purpose of analysis (i.e., radiation modeling) and measurement 
of outgoing flux. <t has been suggested to divide the entire longwave spec- 
tral range into the following subregions (ref. 31): (a) 0.7 to 4 y; (b) 4 

to 8 u ; (c) 8 to 12 u; (d) 9 to 10 u; (e) 12 to 18 u ; and (f) 18 to 50 u . 
Specific reasons for suggesting this spectral subdivision are given in ref- 
erence 33. The Earth radiation budget experiment (ERBE) proposed by NASA 
(ref. 50) consists of two packages designed to provide three spatial resolu- 
tion options with three broad spectral bands as: (a) short wave, 0.2 to 5 u; 
(b) long wave, 5 to 50 u ; and (c) total, 0.2 to 50 y . However, for para- 
metric studies, it is desirable also to extend the long wave range to 200 u 
(i .e. , 5 to 200 u). 

As discussed in the introduction, it is essential to incorporate an 
appropriate model for the cloudy atmosphere in the general radiative trans- 
fer model. For this, information on physical characteristics of clouds is 
essential. Some of the basic information required is: cloud amount, cloud 

height, cloud-top texture, height-width ratios, microphysical properties, 
total water content, liquid/ice water content, and cloud base altitude. All 
these factors contribute to the variability of the cloud radiative proper- 
ties. As such, a simplistic description of clouds, where they are consider- 
ed as opaque/black surfaces, is grossly inadequate. It is essential, there- 
fore, to incorporate into the general radiative transfer schemes appropriate 
cloud models which take into account as many physical variables as possible. 


Basic governing equations for radiative transfer through clouds are availa- 
ble in the literature (refs. 37-40). The upwelling radiance in the presence 
of a cloud layer of thickness Az can be expressed as (ref. 44): 

E(i»),h) = [t£ (Az/T u (4z)]i£ s B(ui,T s )T |u (h,0) 

♦ /o b 8[»,T{z)JdT M (h.z)]} 

+ l" z + 42 B[«,T(z)]dT u (h,z) (2.5) 

b 

Detailed discussions of the above equation and the dependence of cloud 
transmittance too (az) on its liquid-water content and droplet size distri- 
bution is given in reference 44 (a copy of this is attached as Appendix Al). 
In addition, radiation models of finite clouds of different geometrical 
shapes and sizes are becoming increasingly available in the literature 
(refs. 45-49). The existing radiative transfer models should be modified to 
incorporate other cloud parameters in order to investigate the effects of 
their variability on upwalling radiances and fluxes. 

During the International Radiation Symposium held at Colorado State Uni- 
versity in August 1980, it was anphasized that the effect of non-local 
thermodynamic equi 1 ibri in (NLTE) should be considered in the radiative 
transfer formulations for a better interpretation of the Earth radiation 
budget experiment (ERBE) data. Inclusion of this effect was emphasized, es- 
pecially for the fundamental bands of C0 2 , 0 3 , and N 2 0. 

The monochromatic upwelling radiance under the NLTE condition, in the 
presence of a cloud layer of thickness a, may be expressed as (ref. 51): 
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E(«.h) •* [ T ^(A2)/ T <ij (^)]{c s B(w,T s )T a) (h f 0) + E (u Q ) 


+ f? ilnz)dTjh,z)i 

+ ^2 k+ 42 -J.[T(2)]<h B (l>.z) (2.6) 

0 

where J [ T( z) represents the nonequi librium source function and is given 
in terms of the Planck function as 

JjT(z) « { EAu[ T ( 2 ) ] ♦ n R[T(2)}/(1 + n) (2.7) 

where n = n /n r . Equation 2.6 can be ulsed as a diagnostic equation to in- 
vestigate the influence of NLTE on upwelling atmospheric radiation. 

By following the nomenclature adopted by the International Radition Com- 
mission for the Earth Radiation Budget Experiments (ERBE) the expression for 
the upwelling radiance reaching the cloud base from the underlying sur- 
face and atmosphere is given by (refs. 44 and 52) 


MV = 


e B (T )t (z.,0) 
S v v s' v v b* ' 


B v [ T (z)dr v (z b ,z) 


( 2 . 8 ) 


where T v (z b§ 0) is the clear column transmittance with reference to the 
cloud base z^. The radiance reaching the top of the cloud layer at 
z fa + Az is given by 


MV 4z) = M z b )T v T(4z> + !l-\(4;.)]B v ( T c ) 


(2.9) 


where t^(Az) represents the total transmittance through the cloud (i.e., 
it is the product of the transmittance due to cloud and the transmittance of 
the atmosphere in the cloud) and B V (T C ) is the Planck function evaluated 
at the cloud average temperature T . The first term on the right hand side 
of equation (2.9) represents the radiance reaching the cloud top from the 
underlying surface and atmosphere, and the second term is the radiance due 
to the cloud emission. A combination of equations (2.8) and (2.9) results 
in 

M z b Wz) = T J (4Z) e s B v ( V T v' z b’ 0) 

+ * T (az) J 0 b 3“[T(z)]d T (z b , z) 

♦ [1-tJ(Az)]B v (T c ) (2.10) 

The upwelling radiance reaching the top of the atmosphere (denoted by 
h in this study) is expressed as 

L v ( h ) = L v (z b +Az)T v (h, z fa +Az) + J h z B v [T(z)]dT v (h,z) (2.11) 

b 

A combination of equations (2.10) and (2.11) results in 

l v (h) - t'(Az)t v (M)* s 8 v (T s ) 
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+ T 5( A z)/ Q b B v [ T ( z )]dT v (h,2) 


+ T v ( h . z b +Az) f 1 ' T ^ Az} ] B v {T c ) 

+ /i B [T(z)]dx (h,z) 

b v (2.12) 

where t v is the cloud transmittance. It should be noted that equation 
(2.12) is a modified form of equation (2.5). 

c T 

As Az approaches zero, t v (Az) and t v (az) approach unity, and 
equation (2.12) reduces to 

M h > * e sMV T v ( h >°) + l" 0 B v [T(z>] dr u (h,z) (2.13) 

This is an expression for the upwellilng radiance for a plane-parallel clear 
atmosphere and can be expresed also in a direction 6 relative to nadir as 

L (0) = e B (T )t (0) + f 1 B (T ) dT (0) (2.14) 

V v ' s V S VS V J T , v' l' V z ' ' ; 

VR 

In a simplified study, the radiance for an overcast atmosphere is ob- 
tained from equation (2.14) by treating the cloud top as the underlying sur- 
face and considering only that part of the atmosphere which lies above the 
cloud. For partly cloudy conditions, the radiance is obtained as the cloud 
fraction weighted sun of the overcast and clear radiances. 
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The total outgoing radiance is obtained by integrating the spectral 
radiances 

1(e) = / Av L v (0 > dv (2 ' 15) 

By assuming azimuthal symmetry, the spectral outgoing flux can be obtained, 
in principle, by integrating equation (2.15) over the nadir angle 0 from 
0 to ir/2 as 

M = 2ir f 1 ^ 2 L (6) sin 9 cos 0 d 9 (2.16) 

In practice, however, lengthy numerical integration is avoided by using the 
diffusivity approximation (ref. 17). The total outgoing flux M can be ob- 
tained by integrating M v over the frequency range Av as 

M = J A M Mv dv (2.17) 

' Av v 

The anisotropic function is defined as 

R(0) = Tf L(e)/M (2.18) 

The difference between the largest and smallest values of R(0) is a good in- 
dicator of the extent of limb darkening for most atmosphere models. This is 
defined as 

G = R(8 • ) - R (0 ) (2.19) 

v mm ; v max' v ; 

where the largest value of R(0) corresponds to the minimum value of 9 
and vice versa. 
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Figure 2.1. Various components of upwelling atmospheric 
radiation. 
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3. ABSORPTION MODELS 


The greatest problem in computing radiances and fluxes is the integra- 
tion of spectral relations over the frequency range of interest. The 
absorption coefficient (and, hence, the transmittance) is a highly variable 
function of the frequency, and for accurate work it should be evaluated at 
small frequency intervals. Furthermore, within a band which usually con- 
sists of thousands of rotational lines, the absorption coefficient at any 
frequency is made up of contributions from many lines. In principle, there- 
fore, it is possible to calculate the absorption coefficient with very high 
accuracy by summing the contributions of all intervening lines. In prac- 
tice. however, it is a very tedious and time-consuming process. For a wide 
frequency range wit., several bands, each with a large nunber of lines, large 
amounts of computer resources are required. As such, use of simplified 
models for spectral absorption is highly desirable. 

As mentioned in the introduction, considerable efforts have been spent 
in the past in devising simplified models to overcome the problem of numeri- 
cal integration over the complicated line structure of the atmospheric spec- 
trum. A complete review on different absorption models is available in ref- 
erence 4. Hie absorption models proposed for the present study are discuss- 
ed in references 20, 53 and 54. In these references, the validity of the 
QRB model has been investigated by comparing the results of this model with 
the LBL, experimental, and other accurate-model results available in the 
literature. These comparisons snow that the QRB formulation offers an 
accurate and efficient method for calculating the transmi ttances and radi- 
ances in nonhomogeneous nonisothermal systems. Further work, however, is 
needed to establish the validity of this model. It would be desirable to 
compare the results of LBL, QRB, LOWTRAN, K-distribution, and FOSCODE models 
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for selected conditions to establish the validity of the QRB formulation. 
However, because of high computational costs, this is beyond the scope of 
the present study. 
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4. COMPUTATIONAL PROCEDURE AND DATA SOURCE 


The numerical procedure for evaluating the spectral atmospheric trans- 
mittance and the upwelling radiance and the data source used for tfe calcu- 
lations are described briefly in this section. 

In calculating the atmospheric transmittance, the atmosphere is divided 
into a number of layers of equal thickness {in the present case, 1 km). For 
the present study, th a top of the atmosphere was considered first to be 10 
km, which is approximately the top of the troposphere. Selected results 
were also obtained by considering the top of the atmosphere at 20 and 30 
kilometers. The pressure path length is given by the expression 


dU 'j = 5 ij (P j' /P NTP )T NTP /T j )dZ j 


(4.1) 


where Q.^ is the volume mixing -atio of the ith constituent in the jth 
layer, dz. is the thickness of the jth layer, ano and T, are the 

J u J 

pressure and temperature at the center of the jtn layer, respectively . Th e 
transmittance at location z in the atmosnhere is given by 

T(u),r = expf- f u n k..(u>)du ] (4.2) 

j i J J 

Following the procedure for evaluating the atmospheric transmittance, 
upwelling radiance is calculated by dividing the nonhomogeneous atmosphere 
into a nucuer of homogeneous sublayers. The complete numerical procedure 
and computer programs are available in references 53 and 54. 

The line parameters needed for this study {position, strength, line, 
width, etc.) were obtained from McCl atchey et al. (refs. 55 and 56). The 
"McClatchey Tape" is available at the NASA/Langley Research Center. The at- 
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mospheric temperature and pressure profiles were taken from the U.S. Stan- 
dard Atmosphere 1962 {ref. 57). The information on global annual average 
model atmosphere is given in table 4.1. The concentration distributions in 

the atmosphere for H 2 O, CO 2 , N 2 O and O 3 were taken from McClatchey et al. 
(ref. 55). The C0 2 and N 2 0 are assimed to be uniformly mixed in the atmos- 
phere. Rotational and vibrational partition functions, required to account 
for the temperature dependence of the line strengths, were taken from 
McClatchey et al. (ref. 56). 

The categories of other model atmospheres, as obtained from reference 
58, are listed in table 4.2. The information on c ud cover and location is 
given in pertinent chapters. 
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TaMo 4.1 Global annual average model atmosphere 


AlC 

Press 

Temp 

Water Vap 

Ozone 

(km) 

(atm) 

(Kel) 

(PPMV) 

(PPW) 


• .5 

954.61 

284.90 

.695 7E+04 

.2790E-01 

1.5 

845. 59 

278.40 

.5395E+04 

.3078E-01 

2.5 

746.91 

271.91 

.394 9E+04 

.327 7E-01 

3.5 

657.80 

265.41 

.2701E+04 

.3353E-01 

4.5 

577.52 

258.92 

.180 1E+04 

.353 IE-01 

5.5 

505.39 

252.43 

. 1176E+04 

.3891E-01 

6.5 

440.75 

245.94 

.888 7E+03 

.4492E-01 

7.5 

382. 99 

239.46 

.4762E+03 

.5412E-01 

8.5 

331.54 

232.97 

.2692E+03 

.741 8E-01 

9.5 

285.84 

226.49 

. 11 7GE+03 

. 1104E+00 

10.5 

245.40 

220.01 

.5422E+G2 

.170 7E+00 

11.5 

209.84 

216.65 

. 2836E +02 

.2592E+00 

12.5 

179.34 

216.65 

.1534E+02 

.345 IE +00 

13.5 

153.27 

216.65 

. 8613E+01 

.4404E+00 

14.5 

131.00 

216.65 

.5955E+01 

.5726E+00 

15.5 

111.97 

216. 65 

.5940E+01 

. 7369E+00 

16.5 

95.72 

216.65 

.590 3E+01 

.9991E+00 

17.5 

81.82 

216.65 

.5867E+01 

. 1375E+01 

18.5 

69.95 

216.65 

.6291E+01 

.179 6E+01 

19.5 

59.80 

216.65 

. 7358E+01 

.2289E+01 

20.5 

51.13 

217.08 

.901 6E+0 1 

.2793E+01 

21.5 

43.75 

218.08 

. 1151E+02 

.3323E+01 

22.5 

37.46 

219.07 

.147 1E+02 

.3898E+0 1 

23.5 

32.09 

220.06 

. 1867E+02 

• 4391E+01 

24.5 

27.52 

221.06 

.235 5E+02 

.486 7E+01 

25.5 * 

23.62 

222.05 

.2734E+02 

.5371E+01 

26.5 

20.28 

223.04 

.2893E+02 

.590 1E+01 

27.5 

17.43 

224. 03 

.3026E+02 

. 6342E+01 

28.5 

14.99 

225.02 . 

.3154E+02 

.662 8E+0 1 

29.5 

12, 90 

226.01 

. 3276E+02 

. 6673E+01 


Volume Mixing Ratios 
(PPW) 

N20 


C02 

.28 

1.6 

330. 
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Table 4.2 Categories of 106 model atmospheres (ref. 58) 


Cloudiness 

Clearsky cases 

58 


Under cast 

48 

Diurnal 

Daytime cases 


Coverage 

(0700 - 1859) 

52 

(Local Time) 

Nighttime cases 
(1900 - 0656) 

54 

Climate Type 

Maritime 

66 


Continental 
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4. Latitudinal Tropical (0* * 30*) 20 

Distribution 

Mid-lat (30' - 60*) 62 

Sub-arctic polar (60* - 90") 24 

5. Seasonal Spring (3/16 - 6/15) 16 

Coverage 

Summer (6/16 - 9/15) 33 

Autunn/Fall (9/16 - 12/15) 27 

Winter (12/16 - 3/15) 30 
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5. ACCURATE EVALUATION OF LONGWAVE RADIATION IN CLEAR ATMOSPHERE 
In recent studies (refs. 20, 44, 51, and 54), accurate theoretical mod- 
els were developed for evaluating the upwelling atmospheric radiance and 
radiative flux. The existing computer codes were refined to include detail- 
ed information on the LBL and QRB absorption and various surface and atmos- 
pheric parameters. The feasibility of the QRB model for atmospheric studies 
were made to investigate the influence of various surface and atmospheric 
parameters on the upwelling radiances. Some of the important results of 
recent investigations are discussed here very briefly, and the details are 
available in the cited references. 

5.1 Evaluation of Atmospheric Transmittance 
In the longwave range, atmospheric tran quittances were calculated by 
employing the LBL and QRB formulations in selected spectral intervals. The 
calculations were made for clear sky conditions, and contributions of all 
important species were considered in each spectral interval. The results 
were compared with available theoretical and experimental results (ref. 54). 
Specific transmittance results for the spectral range from 2,500 to 2,800 
cm- 1 are shown in figures 5.1 and 5.2. The results calculated by the LBL 
and QRB Mode’s and L0WTRAN program (by considering only the 3.17-p water 
vapor band) are shown in figure 5.1. The results are found to be in good 
agreement. The results presented in figure 5.2 are also for the same spec- 
tral range, but, in this case, contributions of other bands (3.57-g 0 3 , 
3.85-p CH^ , and 4.5-y N 2 0) have been included in calculating the atmospheric 
transmittance. From a comparison of results presented in these figures it 
may be concluded that the QRB results are in good agreement with the LBL and 
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LOWTRAN results and that the transmittance (in this spectral range) is main- 
ly due to the 3.17-u H 2 0 band. Results for other spectral ranges are avail- 
able in reference 54. 

For further study, it is suggested to calculate the atmospheric trans- 
mittance, in different spectral ranges, by employing the QRB, K-distribu- 
tion, LOWTRAN, and FASCOOE models and compare the results with the results 
of LBL formulation. 

5.2 Evaluation of Upwelling Radiance and Radiative Flux 

For the standard atmospheric conditions {table 4.1), the results for 
upwelling radiance and radiative flux were obtained for selected spectral 
ranges by employing the QRB formulation. Ir, each spectral range, contribu- 
tions of all important species were considered. For preliminary study, the 
top of the atmosphere was taken to be at 10 km; but, later results were ob- 
tained also by considering the top of the atmosphere at 20 and 30 km. The 
contribution of the reflected component of solar radiation was includes in 
the calculation of the upwelling radiance for an illustrative case, but the 
contribution of atmospheric radiation reflected from the surface was 
neglected. 

Specific results for the spectral range from 2,500 to 2,800 cm" 1 are 
presented in figure 5.3 for the case in which the contribution of reflected 
solar radiation was included. For a fixed surface emittance, the upwelling 
radiance is seen to increase with increasing surface temperature. This is 
because the surface and atmospheric emissions are relatively higher at high- 
er surface temperatures. For a fixed surface temperature, the radiance is 
seen to increase with decreasing surface emittance. This is because, for 
lower surface emittance, the reflected component of the solar radiation is 
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larger, and this makes the total upwelling radiance relatively higher. 
Similar conclusions were drawn also from the results presented in references 
28 and 54. The results presented in the rest of this section we > 2 obtained 
by neglecting the contributions of the reflected solar radiation. 

For standard atmospheric conditions, the variation in upwelling radi- 
ance and radiative flux at different altitudes is illustrated in figure 5.4 
for the spectral range of 5-50 y where most of the longwave radiative proc- 
esses occur. The results clearly show that the upwelling radiance and flux 
increase with increasing surface temperature and decrease with increasing 
altitude. As pointed out earlier, this, is because the surface and atmos- 
pheric emissions will be relatively higher for higher surface temperatures; 
and the radiative energy will attenuate with increasing atmospheric thick- 
ness. It is further noted that the attenuation above 20 km is relatively 

small, and the top of the atmosphere can be taken as 30 km for most appli- 

cations. 

For different longwave spectral ranges, the variations in upwelling 
radiance and radiative flux with surface temperature are shown in figures 
5.5 and 5.6, respectively. The solid curves represent the results for the 
atmospheric top at 10 km and the broken curves for the atmospheric top at 30 
km. It is seer, that the radiative contributions from different spectral re- 
gions add up to give the highest values for the spectral range 5-50 y; the 
values for other spectral ranges are relatively lower. As noted earlier, 

the results for 30 km are relatively lower than for 10 km. The results for 

the atmospheric window (10.5 - 12.5 y) are the lowest and are identical for 
10 and 3Q kilometers. This is because there are only a few weakly absorb- 
ing species in this range and they do not contribute to the radiative proc- 
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ess beyond certain height. 

The variation of upwelling raaiance and radiative flux with surface 
einittance is shown in figures 5. 7-5. 9. The results presented in figures 5.7 
and 5.8 are for the spectral range 5-50 u, and the results in figure 5.9 are 
for the window region. The results show that for a given surface tempera- 
ture the radiance and flux increase with increasing surface einittance. This 
is true, in general, for all nonreflecting radiating surfaces. For the 
earth-atrnosphere system, however, the trend exhibited in figures 5.7 - 5.9 
is true only if the contribution of reflected solar radiation is neglected. 

From the results presented in this- section, it may be concluded that 
the QRB formulation is quite suitable for most atmospheric applications. 

The procedure for calculating the upwelling radiance and radiative flux has 
been developed by using the QRB model and illustrative results have been ob- 
tained for the standard atmospheric; conditions. The results demonstrate 
that for a fixed surface emittance the upwelling radiance and flux increase 
with increasing surface temperature and decrease with increasing height. 

The radiative attenuation above 20 km is small and the top of the atmosphere 
can be considered as 30 km for most applications. The radiative contribu- 
tions from different spectral ranges add up to give the total upwelling 
radiance. The contribution from the window region is quite small and the 
spectral region 5-20 contributes considerably to the entire longwave radia- 
tion. 
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Figure 5.1 Comparison of atmospheric transmittance in the 

spectral range from 2500 to 2800 cm' 1 considering 
water vapor only. 
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Figure 5.2 Comparison of atmospheric transmittance in the 
spectral range from 2500 to 2800 cm" 1 • 
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Figure 5.5 Variation in upweliing radiative with surface 
te: ~>erature for different spectral ranges. 
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6. SENSITIVITY ANALYSIS OF UPWELLING RADIANCE IN PRESENCE OF CLOUDS 

The existing radiative transfer programs were modified to investigate 
the effects of clouds on the upwelling radiance. In order to study the sen- 
sitivity of upwelling radiance to variations in cloud height, cloud liquid 
water content, and cloud thickness, calculations were made in three differ- 
ent spectral ranges. They are the 5- to 10-um, 10.5- to 12.5-ym (window), 
and 10- to 2Q-yrn regions. These spectral regions include most of the strong 
absorption and emission bands of 1^0, C0 2 , 0 3 , and N 2 0. The QRB formulation 
was employed and results were obtained for the standard atmospheric condi- 
tions. The key results are discussed completely in reference 44, and a copy 
of this is attached as Appendix Al; the entire results are tabulated in 
Appendix B2. 

As an example of the sensitivity analysis, the variation of radiance 
with cloud height is illustrated in figure 6.1 for different spectral range 
and surface emissivity. It is seen that the upwelling radiance decreases 
with increasing cloud height for all three spectral ranges. This is because 
the effective cloud top temperature at which the absorbed radiance is re- 
emitted is lower for higher level clouds. It is noted that the maximum im- 
pact of the cloud layer is in the window region. Here, the reduction in up- 
welling radiance (from the clearsky value) by a cloud at z = 11 km is ^bout 
49 percent as compared to about 43 percent in the 5- to 10-u range and 37 
percent in the 10- to 20-um range. 

Other results of this study presented in Appendix Al indicate that the 
difference between the clear and cloudy sky radiances is nearly 36 percent 
for the cloud base at 5 km and about 43 percent for the cloud base at 10 km; 
this difference increases with increasing surface emissivity and surface 
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temperature. Further sensitivity studies, however, are needed to study the 
formation of different types of clouds ar. i tneir interaction with the atmos- 
phere, study the influence of cloud height variations on the net upwelling 
flux, evaluate the total longwave flux and correlate this with the window 
flux by varying the cloud droplet-size distribution, and evaluate the vari- 
ation of total and window region flux with the cloud liquid water content 
and cloud top temperature. 
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7. EVALUATION OF ANISOTROPIC FUNCTIONS IN THE LON&IAVE REG’.O.. 

7.1 Introduction 

A study of the anisotropic function in the atmosphere has ass lined 
greater importance recently because of its applicability for the Earth Radi- 
ation 8udget Experiment (ERBE). This experiment is designed to measure the 
radiation budget of the Earth- atmosphere system at the top of the atmos- 
phere. The scanning radiometer which is an important component of the 
instrument complement measures the radiance emanating at the top in a given 
direction. These directional radiances can be converted to the total out- 
going flux density only if the anisotropic functions for the atmosphere are 
known. 

In this study, anisotropic functions have been calculate^ for several 
model atmospheres. The effects of variability of various meteorological 
parameters on the anisotropic function are investigated. During the initial 
phases of this study, the spectral range considered for the longwave radia- 
tive transfer was 5-50u. The results were presented at a national confer- 
ence (AIAA Paper 83-0161) and the entire material is included in this report 
as Appendix A2. Some specific results discussed in Appendix A2 are tabula- 
ted in Appendix 83. Later, the computer code was modified to cover the 
longwave spectral range 5-200y. Key results of this study are presented 
here and are compared with the results presented in Appendix A2. 

7.2 Physical Conditions and Computational Procedure 

The physical conditions, data source, and computation procedure are es- 
entially the same as given in Appendix A 2 ; these are discussed here 
briefly. 
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The top of the atmosphere is considered to be at 30km. For calculation 
of the transmittance and radiance, the atmosphere is divided into 15 
layers; the first ten layers are one-kilometer thick. The low altitude 
cloud is assuned to be at 2 km, middle-level cloud at *6 km, and high-level 
cloud at 10 km. The thickness of the cloud layers is assuned to be one 
kilometer 

Several cl imatological-average model atmospheres are used to establish 
the latitudinal and seasonal variability of the anisotropic functions. 

These are tropical, mid-latitude simmer/ winter, sub-arctic simmer/ winter in 
addition to the U.S. Standard Atmosphere (USSA) which represents the mean- 
annual mid-latitude conditions. Pressure, temperature, and ozone profiles 
for the above models are taken from McClatchey et al. {refs. 55 and 56). A 
climatological mean value for the surface relative hunid :y (RH) of 75% 

(ref. 59) is used for all of the above models and verti il distribution of 
water vapor is computed using the power law (ref. 60) 

w z = * z ( p z /p s )X t 7 - 1 ) 

where A is related to the water-vapor scale height and its average value 
is taken to be 3. The concentration of water vapor is calculated by 

C = A exp(18. 9766-14. 9595A-2. 438»\ 2 ) (7.2) 

where A = 273. 15/T^. Upon dividing by the density and multiplying by the 
ratio of air molecular weight to the water vapor molecular weight, the con- 
centration of water vapor in ppmV is obtained. Data files are created for 
different surface relative humidities. 

The computer code developed for this study has the following capabil- 
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ities: 


1. The upwelling radiance and radiative flux in the longwave range from 5 
to 200 u {2000 to 50 cm -1 ) can be calculated; the program can be used 
also in any specified spectral range between 5-20CT u. 

2. The participating species considered for the calculation of the atmos- 
pheric transmittance are Ch 4 , C0 2 , H 2 0, N 2 0, and 0 3 . The effect of any 
particular gas (i.e., the radiative contribution of a specific molecu- 
lar gas) can be investigated. 

3. The atmosphere (with top at 30 km) is divided into 15 sublayers; the 
first ten kilometers is divided into ten equal layers. The cloud base 
can be considered anywhere below 10 km. It is assumed that most clouds 
lie below 10 km. 

4. The parameters that can be changed easily in the program are surface 
temperature, surface emissivity (emittance), cloud height and cloud 
emi ssivity. 

5. Radiances and fluxes can be calculated for different cloud-cover frac- 
tions and zenith angles. 

6. Anisotropic functions can be calculated for different atmospheric mod- 
els. 

For a given atmospheric model, the computer program is capable of com- 
puting all results (including all information given in steps one through 
six) in one step. This reduces the computer time considerably as compared 
to the individual runs. For parametric and sensitivity studies, therefore, 
the present code is very economical. The information on the computer 
program is available in Appendices Cl and C2; Appendix Cl provides the 
listing of the program. Comment cards are inserted in the listing to ex- 
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plain each stage of the program. 


7.3 Results and Discussion 

The sensitivity of the anisotropic functions to meteorlogical variables 
like cloud cover, cloud height, surface relative humidity, surface emittance 
and temperature is examined for the U.S. Standard, Tropical, and Sub-Arctic 
Winter atmospheres. 

Tlie results of clearsky radiance and flux are tabulated in table 7.1 
for the U.S. Standard Atmosphere for different longwave spectral ranges. 
These results are similar to those presented in chapter 5; but now the spec- 
tral range has been extended to cover 5-200y. The results show that the ra- 
diance and flux for 5-200 y range are only about three-and-one-half percent 
higher than the values for 5-50 y range. This is not a significant in- 
crease, but this information may be essential for some specific applica- 
tions. 

For the standard atmosphere and spectral range 5-200u, figure 7.1 shows 
the upwelling radiance as a function of the cloud cover fraction for differ- 
values of the zenith angle. The cloud top is assumed to be at 5 km and the 
cloud emissivity is taken to be unity. As would be expected, the results 
clearly show that the upwelling radiance decreases with increasing cloud 
cover and zenith angle. The radiance values for 0=0° and 15* are seen to 
be quite close, but considerable decrease in radiance is noted for 0-values 
larger than 45°. As explained in Chapter 2, the difference between the max- 
imum and minimum values of radiance represents the extent of limb darkening 
for most atmospheric models. 

The significance of the anisotropic function and importance of the 
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Table 7.1 Clearsky radiance and flux for U.S. standard (mid-lat. 

average] atmosphere. Atm Top = 30 km, E = 1.0, T = 288.15 K 

s s 

surface relative tumidity (RH) = 7 5%, 6=0. 


Spectral Range, 

V 

Clear Atm. 
Radiance, W/rn^sr 

Clear Atm. Flux, 
W/m 2 

5-10 

16.63 

45.58 

10 - 20 

42.16 

124,94 

5 - 20 

58.79 

170.52 

10.5 - 12.5 

14.79 

45.58 

5 - 50 

78.59 

229.10 

5 - 200 

81.41 

237,85 









limb-darkening work are illustrated in figures 7.2 and 7.3. It is noted 
that the radiance emerging from the top of the atmosphere can be highly ani- 
sotropic (non Lambertian). For a given Earth- atmosphere system, if the radi- 
ation leaving the atmosphere is isotropic (Lambertian), the satellite in- 
strunent will receive the same radiation for all nadir viewing angle. How- 
ever, because of the directional dependence of the upwelling radiation, dif- 
ferent radiative energy is received by the instrument at different viewing 
angle. The extent of the limb darkening is expressed by the value of G 
which is the difference between the maximum and ininmum values of the aniso- 
tropic function for a given atmospheric model. The results for anisotropic 
functions are tabulated in Appendices 83 and 84 for the spectral ranges 
5-50 v and 5-200 u, respectively. The results for the spectral range 5-50 y 
are discussed in Appendix A2 and some specific results for the spectral 
range 5-200 u are discussed in this section. 

The latitudinal variability of the anisotropic functions for the clima- 
tological-average model atmospheres is shown in figure 7.4; the values of G 
for different models are given within the parentheses. In comparison to the 
results for the mid-latitude average model, the values of the anisotropic 
function are found to be ' igher for the tropical atmosphere and lower for 
the sub-arctic winter atmosphere for lower nadir viewing angles; however, 
the reverse trend is observed for vievring angles greater than 50°. The val- 
ue of G for the mid- latitude average clevr atmosphere was found to be 0. 2386 
(instead of 0.2444 for the spectral range 5-50u) and is used as a reference 
value in the following discussions. The values of G for the tropical and 
sub-arctic winter models were found to be 0.3073 and 0.1622, respectively. 
Similar latitudinal variability in anisotropic functions was observed when 
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radio-sor.ae-measured model atmospheres were used instead of the climatolog- 
ical-average models (see Appendix A2). 

The sensitivity of the isotropic functions to changes in the values of 
various meteorological parameters was examined in detail for the U. S. stan- 
dard (mid-latitude average), tropical and sub-arctic winter atmospheres. 

For all case considered, the top of the atmosphere was taken to be 30 km. 

The entire results are presented in Appendix B4, selected results are given 
in tables 7.2 and 7.3, and some specific results are illustrated in figures 
7. 5-7.8. 

The effects of cloud height on the anisotropic functions were examined 
for the overcast (100% cloud cover) conditions and the results for the stan- 
dard and tropical atmospheres are shown in figures 7.5a and 7.5b, respect- 
ively. It is seen that R(0) decreases with 0 sharply for the lower 
cloud heights than for the cloud at z c = 10 km; the variation is steeper 
for the tropical atmosphere than the standard atmosphere. The value of G 
varies from 0.2013 (z c = 2 km) to 0.0186 (z c = 10 km) for the standard 
atmosphere and from 0.2646 (z * 2 km' to 0.1050 (z = 10 km) for the 

V w 

tropical atmosphere. The results clearly indicate that the anisotropic 
atmosphere functions are quite sensitive to the location of clouds in the 
atmosphere. 

The sensitivity of the anisotropic function to cloud covers is illus- 
trated in figures 7.6a and 7.6b for the standard and tropical atmospheres, 
respectively. The results are obtained for a cloud height of 6 km and the 
cloud emissivity is taken to be unity. It is noted that the value of G 
varies from 0.2386 (clear sky) to 0.1265 (overcast sky) for the standard at- 
mosphere and from 0.3073 (clear sky) to 0.1958 (overcast sky) for the trop- 
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Table 7.2 Sensitivity of anisotropic functions to various meteorological 
parameters for a mid-latitude average atmosphere. Atm Top = 

30 km. 


No. 

Meteorological Parameter 

G-values 

5-50 y 

5-200 y 

1 . 

Surface relative humidity, % 




50 

0.238 



75 

0.244 

Wmmm 


ICO 

0.249 

HU 

2. 

Water vapor scale-height 




parameter, X 




2 

0.258 



3 

0.244 



4 

0.233 


3. 

Surface emittance, t 




0.8 

0.201 

0.1966 


0.9 

0.223 

o no. 


1.0 

0.244 

C 

4. 

Cloud height (overcast), km 




2 

0.205 

0.2013 


6 

0.128 

0.1265 


10 

0.014 

0.0186 

5. 

Cloud cover (z c = 6 km), % 




0 

0.244 

0.2386 


50 

0.195 

0.1906 


100 

0.128 

0.1265 

6. 

High-cloud emissivity, e c 




0.5 

0.162 

0.1586 


1.0 

0.014 

0.0186 














ical atmosphere. The results show that the anisotropic functions are also 
quite sensitive to the changes in fractional cloud cover. 

The sensitivity of the anisotropic function to the high-cloud emissiv- 
ity is shown in figures 7.7a and 7.7 b; the results are’ obtained for an over- 
cast sky. The value of G is seen to vary from 0.1586 (e = 0.5) to 0.0186 

(e - 1.0) for the standard atmosphere and from 0.2318 (e = 0.5) to 0.1050 
c c 

(e = 1.0) for the tropical atmosphere. The anisotropic function is seen to 
be quite sensitive also to the high-cloud emissivity. 

The strong dependence of G on cloud height is a combination of two 
effects. Firstly, the lower temperature of the cloud top reduces the 
temperature difference between the underlaying surface (cloud top) and top 
layers of the atmosphere. Secondly, there is much less water vapor above 
the cloud top. Variation of G with fractional cloud cover is simply a 
combination of the effects for clear and overcast cases. The large increase 
in G as the high-cloud emissivity decreases from 1.0 to 0.5 can be attri- 
buted to the fact that for the partly transmitting cloud, considerable part 
of the '■adiation emanating at the top originated at the surface. It should 
be noted that the values of G for 50% cloud cover and = 1.0 is equal 
to chat for 100% cloud cover and = 0.5; this demonstrates that cloud 

emissivity and fractional cloud cover are equivalent parameters (i.e., the 
effective cloud cover is e c x actual cloud cover). 

The sensitivity of the anisotropic function to the surface erittance 
(emissivity) is shown in figures 7.8a and 7.8b; the results are obtained for 
the clear-sky conditions. The value of G is found to vary from 0.1966 
(e * 0.8) to 0.2386 (e $ = 1.0) for the standard atmosphere and from 0.2590 
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(e = 0.8) to 0.3073 (e $ =1.0) for the tropical atmosphere. It is noted 

that the anisotropic function is also quite sensitive to the surface emis- 
sivity; an increase in surface emittance of 10% results in an increase of G 
value of about 11% (see table 7.2). The sensitivity of the isotropic func- 
tion to the variation of surface temperature was examined for different mod- 
els (see Appendices A2, B3, and B4), A 5 K change in the surface tempera- 
ture (without changilng the temperature profile) was found to cause more 
than 6% change in the G value. It should be noted that the changes in 
surface emittance and temperature effectively amount to a change in the 
lapse rate and this, in turn, causes the change in G values. When atmos- 
pheric temperature profile was changed along with the surface temperature, 
insignificant changes in G were observed (Appendix B3). 

The sensitivity of the isotropic function to the variatiun of surface 
relative hur.idity (at fixed surface temperature) water vapor scale-height 
parameter A, and CO 2 concentration was examined for different models (see 
table 7.2 and Appendices A2, B3, and B4). It was found that G increased 
sharply from 0.091 for a dry atmosphere (no water vapor) to 0.221 for 5% 
relative tumidity (/fcpendices A2 and B3), but increased very slowly there- 
after. The variation of the scale-height parameter was found to have a 
small effect on G (table 7.2). The variation of carbon dioxide concentra- 
tration between zero and twice the standard amount ( i . e. , 660 ppmV) has very 
small effect on G (Appendix B3). 

7.4 Conclusions 

The existing computer code for the 5-50 u longwave range was modified 
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to cover the spectral range 5-200 p. Other modifications were made to 
achieve higher efficiency and less computational costs. The revised code 
was used to investigate the effects of various meteorological parameters on 
the upwelling radiation and anisotropic function. The study shows that in- 
clusion of the spectral range from 50 u to 200 m amounts to about 3% change 
in final results in most cases. This may not appear to be a significant 
change- but for some specific applications this information may be of vital 
importance. Other conclusions on "Evaluation of Anisotropic Functions in 
the Longwave Region" are essentially the same as given under "Concluding Re- 
marks" of Appendix A2. 
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Figure 7.3 Significance of the anisotropic function and 
definition of G. 
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Figure 7.4 Latitudinal variability of the anisotropic functions 
for the cl imatological-avorage model atmospheres. 
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Figure 7.5a Anisotropic functions for selected values of cloud-top 
height, U.S. standard atmosphere. 
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Figure 7.Eb Anisotropic functions for selected values of cloud-too 
neight, tropical atmosphere. ' ’ 
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Figure 7.6b Anisotropic functions for different values of 
fractional-cloud cover, tropical atmosphere. 
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8. CONCLUDING REMARKS 


Radiative transfer models and computer codes have been developed by em- 
ploying the line-by-line and quasi-random band formulations to determine the 
gas emissivities, upwelling atmospheric radiance and radiative flux in the 
entire longwave spectral range (5-200 u). The program is quite versatile 
and can be used to investigate the influence of various parameters and non- 
equilibriim radiation on the net radiative exchange. 

In this study, the program was used to evaluate the radiative flux in 
clear atmosphere, provide sensitivity analysis of upwelling radiance in 
presence of clouds, and determine the .effects of various climatological par- 
ameters on the upwelling radiation and anisotropic function. The results 
show that the QRB formulation is quite suitable for most atmospheric appli- 
cations, the top of the atmosphere can be considered at 30 km, and the radi- 
ative contribution from the spectral range 50-200 u anounts to about three 
percent. For a given surface and atmospheric conditions, the upwelling 
radiance and anisotropic function are found to be very sensitive to the 
variation in cloud parameters (liquid-water content, thickness, height, 
cover, and emissivity). Other studies, however, are needed to further in- 
vestigate the specific influence of various cloud parameters. It is estab- 
lished that the 1 imb-darkeni ng in the atmosphere is caused primarily by the 
presence of water vapor and reaches saturation for very low values of water 
vapor burden. 
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APPENDICES 


APPENDIX B1 

CLEARSKY UPWELLING RADIANCE AND RADIATIVE FLUX 
TABLES Bl.l - B1.5 
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Table B1.2 Clearsky upwelllng radiance and radiative flux for spectral range 5-50 p. 
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Table B1.3 Clearsky upwelllng radiance and radiative flux for spectral range 10-20 u 
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Table Bl.o Clearsky up welling radiance and radiative flux for spectral range 10.5-12.5 u 
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APPENDIX B2 

UFWELLING RADIANCE IN PRESENCE OF CLOUDS 
TABLES B2.1 - B2.4 


Table B2.1 Upwelllng radiance In presence of clouds [(W/cn^-sr) x 10” 5 ] 

for V = 280 K and spectral range = 5-10 u (1,000 - 2.000 cm -1 ) 
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Clearsky Liquid Water Content (LWC), gra/in 2 


OF POOR QOhU«t 


ja 

pm e 


H O O Cft CO N 

Nooch^o^ 

N<»hN^<h 
w UH6 iO iv 00 

io n c n n o 

IT CO IS s (30 O) 

CO G1 o c\i ro 

HH^(SjCd(S1 

OO O' O ^ 

— i — ' CM CM CM CM 

Wtofsco^o 

^ *■( »^ M ^ CM 

o ma) 0 U)« 

anpsrtioot 

w-irstM 

rsO«co<-iifl 

m .— « rs. cm co 
CO CM in O' CM lO 

^ cm co in in 
CM CM CSJ CM C\J CM 

h fo ^ in r« aj 

CM CM CM CM CM CM 

oo o *-h cm «a* m 

CM CM CM CM CM 

ao o 

o> P«- «*• cm o 

mN^NaiCM 
CO 'O CO O ISO) 

OO 10 CO ^ ^ ID 
ON^CMOIM) 

IAiO00ON«r 
cm cm cm m eo <*i 

m fs <7\ — • cm ■«r 

CM CM CM OO CO O") 

fO^lOCOOl-H 
CM CM CM CM CM C"> 

<y» m oo cm f*s 

»>»NtN 
• ••«•• 
O CNJ ^ N 0\ ^ 
m co ro ^ 

N VOtO 
^ co m r*. as 
«•#••• 
—« i 1 ) to rs* cn cm 
fO or ro ro ro ^ 

in ^ ro co cm csi 
cnmrsoiHPO 
#••••• 
flOOCM^-rso, 
cm co co co ro co 

o Ln o tn h 
cm a> co vo m 

oo lo m ^ cm 

00 VO ^ CM O 00 

»-♦ oi co n m 
— 1 oo vo *#■ !M o 

rs o cm m co 
oo ^ *r m 

NOMVOOIH 
po xf tj ^ m 

insOfiioo) 

n voco wco vo 

CO ^ O N CO Ol 

CM CO — i Ln Oi ■M" 
♦ oun Oioh 

o' ^ ao <*■> oo cm 

lOCMNdOO* 

ui wtooofo 
m m 10 10 

coo'miso^ 
*#■ in m lo lO 

O’} IS o ^ Is ^ 

^«u)imoio 

^ in vo r** ao oi 
n co Ci io »h rs. 

(O IS lO 00 O' o 
CM T' CM IS CM CO 

HMCM^^in 

VO 10 f-H 10 t-H 

OflU)O^Cl 
lO O >C N Ms 

Is -no Olfl Ol 
WCOcOCsisN 

«t O' CO CO CM IS 
ininiotoNis 

•-* 00 ** O m CM 
'C^COCSJOO 

WOOh^isO 

OSlrtCSJOlN 

cti on in co *-< *• 

* CM 0"O <■ f 

o io m oo 't w 
rs rv 00 00 O' 

CO CM 00 CO O' 

n ts no oo O' <Ji 

ao 4 1 O' m fh is 
lOlsNQO O' O' 

i^QUlO Wh 
'MNH toOlO 

OOCOIOCOhO 

in oo — * ^ is. o 

(O cm ^ c?> oo r-s 
n cn cm m 

oOinnoO'n 

oo o> O *-t CM 

*■4 H ^-4 

cow^onu) 

mOlOHHCVi 

iH H H *H 

m no O rs in cm 
ao O' O a f-< cm 

f~F .H fH i ( 

awi’i^ioh. 

lOO'flOM'O 

cvj m m io 

^hi n in d h 

M nl iniO Is 

eo m rs O' f-i 

KN^OlOUO^- 
O •— l CM ro ^ LT> 
^ H H H H H 

aors-^ninTf*^* 
o ^ cm cn ^ ^n 

^ H H ^ 

co is io m ^ 
OnNcof in 

^4 tH H pH 

LO O PsCO ^ O 

in co rs oo o> o 

m io is • go ci o 

o o o o o — « 

dodooH 

O O O O O h 

o 

o 

.*-4 

o 

CM 

o 

CM 

O 

CO 

o 

ro 




















ORIGINAL FASt l* 
OF POOR QUALITY 


O a 

— . « o 

*- CM ■«- 

I/) I _l 

I O 


OO^^hCM 

^ONCOOiOh 

(Njcsi 

H ^ H fH 

H ^ H ^ H ^ 

m o vo n o vo 
tn to VC r. 00 CO 

to vo oo O' o 
• * • • ■ • 
io id to id vo 
0 9l O' 01 O' 0) 

95.879 

95.988 

96.097 

96.206 

96.315 

96.423 

O^CON ^Oio 
^ P*"* 0 H fO 

• ••••• 
CM CM CM CM CO CO 
H ^ H H ^ 

H H »■*< ^ ^ 

msMtoO^ 
u» p tn <*) ey 
nin rv (ji h to 
• ••••• 
NNh-NCOCO 
01 O' O' Oi 01 01 

^ CM .-1 O' 00 CO 
Ot O' O' 00 00 00 

VO 00 O <NI •M- VO 

lOlONNNN 
01 01 01 01 01 01 

r**> CM vO ■— • *n O 
ao cm m cn cm lo 
• ••••• 

^ ^ H ^ H i-4 
H ^ H ^ H fH 

Hinoconoo 

hiohnmi/> 

00 «— • VO 00 • • 

. • . . o o 

COOiOIOtOO 

01 W Ot Ot H H 

0 i «-h cm ^ vo 
HT O 40 CM 00 
i-4 10 01 CM IQ 01 
■ ••••• 
00 00 00 01 01 01 
01 Ol 01 0 01 0 

TfPv O to lOOl 
U)'<C0'tO'O 

NHiootnot 

^Hfs^OtO 

O' 4f O VO CM 00 

oo (A cm ao in 

^ONn00 9iOI 

H »H H «H pH ^4 
pH fH pH fH pH pH 

HNNIO^vf 

c o o o o o 

^H »H fH ^ ^H v— 4 

O HN N CO p 
O o O O O O 

rH r4 ^4 .4 i—4 

NIOh«OOIA 
^iflrscoo^ 
*•!••• 
f-t N n ♦ «j 
CM CM CM CM CM CM 
*""* (H H4 *H 

HCO VOCOHOO 
VO 00 O H 
• « • • • • 
VONCOOIHN 
O O O O 
H H H i — 4 H »■ "4 

CM CM CM CM CM CM 
CO O CM vf vO CO 
...... 

vf> 1 — go Cfl o •— 1 

0 0 0 0.-4.-4 
. 4 M ^ 

CT«aor«.cOLn^f 
cm n ^ in vo i — 

rti>.HtflOt(*> 

HCV^tflVOOO 

p Ot CO to VO O’ 
SI OtHlOION 

o CM ^ VO 00 o 

rt n n fl M vf 

pH fH ^4 H ^H pH 

ION Ot HM 1ft 
HHHOJNN 
H W H I— 4 H i-C 

^'OO'Hfoui 
pH «H phCM CM M 
H ^H pH *H ^H pH 

VO N 00 0 o o 
^NOCONin 

O O O O O O 
{MHOOlCOS 

1 — vO VO co CM O 

^H p“H pH pH ^H pH 

«0 ^rsi ^in 

io to vn vo vo 

H ^H ^ pH pH 

HlflWCMVOO 
<n n n vn 

^ pH ^H pH p-H pH 

HIOOKOI^H 
(O <0 (O vt «t Sj 

*— • #H i—4 rH *—4 — 4 

rtNHlflOtN 
O' 00 CO P* VO VO 

» Ol O O h (VI 
*t to r*> 00 O' O 

O C*^ CO 0 VO CM 

o cm m o co 

to cm ov vo no 

psCOCO® OH 
CM 

ON^HOVO 
VO VO P* 00 00 Ov 

^4 H r-4 H .-4 i— 1 

-iooioNON 

VO VO 1 — CO O' O' 

pH fH pH f— 4 pH H 

PHNCM^O^ 
CO CM o> m CM CO 
vO • • • • • 

• CM ^ o CNJ 

0 ^ in vo co 0 

CM CM CM CM CM CM 

MJiPOiP O 
00 r-^ vo vo vo 

...... 

n vo O' cm vo oo 

•— t CM oo to CO 1 — 
CM CM CM CM CM CM 

215.36 

228.59 

241.83 

255.06 

268.30 

281.53 

310.49 

332.89 

355.29 

377.69 

400.09 

422.09 

296.34 

319.31 

342.28 

365.25 

388.23 

411.20 

^ fH Q0 iil pH CO 
coor^cvjiN. h 
• ••••• 
0 C0\0 o CON 
O'CVI^NO'H 

cm no co <v-> ^ 

vo 'O co a o 
...... 

o o o o o 

irMONODOlO 

• •■•if 

0 0 0 0 0—4 

in vo p*. oo o' O 
• ••••• 
O 3 O O O 

in 

VO 

m 

o 

pH 

o 

CM 

o 

CO 


[( 



















Table 82.2 (continued) 
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Table B2.3 Upwelllng radiance In presence of clouds [(W/cm 2 -sr) x 10 -5 
for T 5 * 280 K and spectral range - 5-10 g (1,000-2,000 cm" 
















Table B2.4 Upwelling radiance In presence of clouds [(W/cn^-sr) x 10“ 5 ] 

for T = 280 K and spectral range = 10.5 - 12.5 p (800-950 cm -1 ) 
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Table B2.4 (continued) 
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APPENDIX B3 


ANISOTROPIC FUNCTIONS FOR SPECTRAL 
RANGE 5-50 - Tables B3.1 - B3.12 


Table B3.1 Variation of anisotropic function with cloud cover fraction 

for standard atmosphere; Atm Top * 30 km, spectral range = 5-50 
V, T $ = 288,15 K, e s = 1.0, \ = 3, e r = 1.0, Z r = 6 km, RH 


Zenith 
Angle, 0 

Cloud Cover, % 

0 

50 

100 

0 

1.079 

1.062 

1.040 

15 

1.073 

1.058 

1.037 

30 

1.056 

1.044 

1.029 

45 

1.024 

1.019 

1.013 

60 

0.972 

0.978 

0.985 

70 

0.919 

0.935 

0.957 

80 

0.835 

0.867 

0.912 

G 

0.244 

0.195 

0.128 


101 








Table B3.2 Variation of anisotropic function with cloud height for 
standard atmosphere; Atm Top = 30 km, spectral range 
= 5-50 u, T $ = 288.15 K, e s = 1.0, X = 3, e c = 1.0, RH = 75% 


Zenith 



Cloud Height (Z c ), km 








Angle, 

9 

50% Cloud Cover 

100% Cloud Cover 

2 

6 

10 

2 

6 

10 

0 

1.073 

1.062 

1.092 

1.067 

1.040 

1.005 

15 

1.068 

1.058 

1.049 

1.063 

1.040 

1.005 

30 

1.052 

1.044 

1.037 

1.048 

1.029 

1.004 

45 

1.023 

1.019 

1.016 

1.021 

1.013 

1.002 

60 

0.974 

0.978 

0.982 

0.976 

0.985 

0.998 

70 

0.925 

0.935 

0.946 

0.931 

0.957 

0.995 

80 

0.848 

0.867 

0.890 

0.862 

0,912 

0.991 

G 

0.225 

0.195 

0.162 

0.205 

0.128 

0.014 
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Table B3.3 Variation of anisotropic function with high-cloud emissivity 
for standard atmosphere; Atm Top = 30 km, spectral range 
= 5-50 y, T $ = 288.15 K, e $ = 1.0, X = 3, RH = 75% 


Zenith 

Angle, 

0 

High Cloud (Z c * 10 km) 

Emissivity, e c 

50% Cloud Cover 

100% Cloud Cover 


0.5 

1.0 

0.5 

1.0 

0 

— 

1.052 

— 

1.005 

15 


1.052 


1.005 

30 


1.037 


1.004 

45 

1.021 

1.016 

1.016 

1.002 

60 

0.976 

0.982 

0.982 

0.998 

70 

0.931 

0.946 • 

0.946 

0.995 

80 

0.859 

0.890 

0.890 

0.991 

G 

0.208 

0.162 

0.162 

0.014 
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Table B3.4a Variation of anisotropic function with surface relative 

hunidity for standard atmosphere; Atm Top = 30 km, spectral 
ranges. 5-50 u, T c * 288.15 K, e s = 1.0, X = 3, e r = 1.0, 


Zenith 


Surface Relative Hunidity (RH), % 


Angle, 

e 


Clear Sky 


Overcast sky 


75 100 25 50 75 100 


075 1.077 1.079 1.0 

069 1.072 1.073 1.0 

053 l .055 1.056 l.fl 

023 1.024 1.024 1.0 

974 0.973 0.972 0.9 

923 0.921 0.919 0.9 

0.843 0.339 0.835 0.831 


0 

0 

0 

0 

9 

9 


0.914 


.039 
.036 
,.028 
.012 1.013 1.013 
1.986 0.905 0.985 
1,957 0.957 0.956 
0.912 0.912 0.912 


1 . 

040 

1 . 

037 

1 . 

029 

1 . 

013 


905 


957 


0.238 0.244 0.249 0.122 0.127 0.128 0.129 































Table B3.4b Variation of anisotropic function with surface relative 
humidity (low range) for standard atmosphere; Atm Top 
= 30 km, spectral range = 5-50u, T $ * 288.15 K, 

£ g = 1.0, \ = 3, clear sky 


Zenith 

Angle, 

e 


Surface 

Relative Humidity (RH), % 


0 

5 

10 

15 

20 

25 

0 

1.028 

1.068 

1.071 

1 073 

1.074 

1 .075 

15 

1.026 

1.063 

1.066 

1 068 

1.069 

1.069 

30 

1.020 

1. 37 

1.051 

1 052 

1.053 

1.053 

45 

1.009 

1.022 

1.022 

1.023 

1.023 

1.023 

60 

0 990 

0.375 

0.9/4 

0.974 

0.974 

0.974 

70 

0.970 

0,926 

0.924 

0.924 

0.923 

0.923 

80 

0.937 

0.847 

0.845 

0.844 

0.843 

0.843 

G 

0.091 

0.221 

0.226 

0.229 

0.231 

0.232 



I 


i able B3.5 Variation of anisotropic function with water vapor scale-height 
parameter for standard atmosphere; Atm Top = 30 km, spectral 
range = 5-50 u, T $ = 288.15 K, c $ = l.u, l Q = 6 km, RH = 75% 


Zenith 

Angle, 

0 


Water Vapor Scale-Height Parameter, > 


Clear Sky 



2 

3 

0 

1.084 

1.079 

15 

1.078 

1.073 

30 

1.060 

1.056 

45 

1.026 

1.024 

60 

0.971 

0.972 

70 

0.915 

0.919 

80 

0.826 

0.835 

G 

0.258 

0.244 


1.074 

1.069 

1.053 

1.023 

0.974 

0.923 

0.841 


■ 


..044 

.041 

.031 

.014 

0.985 

0.956 

0.914 


0.233 0.130 


Overcast Sky 


1.040 

1.037 

1.029 

1.013 

0.985 

0.957 

0.912 


0.128 
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Table B3.6 


Variation of anisotropic 
standard atmosphere; Aim 


= 5-50 u, T s = 288.15, \ 


function with surface 
Top = 30 km, spectral 
= 3, e c = 1.0, Z c = 6 


emissivity for 
range 

km, RH = 75%. 


Zenith 

Angle, 


Surface Emissivity, e s 

- 



Clear Sky 


50% Cloud Cover 

0 

0.8 

0.9 

1.0 

0.8 

0.9 

1.0 

0 

1.066 

1.073 

1.079 

1.055 

nn 

1.062 

15 


1.068 

1.073 

1.051 

| mm 

1.058 

30 

1.047 

1.052 

1.056 

1.039 

WSsM 

1.044 

45 

1.021 

1.023 

1.024 

1.017 

EEs 

1.019 

60 

0.977 

0.975 

0.972 

0.981 


0.978 

70 

0.933 

0.926 

0.919 

0.944 


0.935 

80 

0.865 

0.850 

0.835 

0.886 

KbB 

0.867 

G 

0.201 

0.223 

0.244 

0.169 

0.182 

0.195 
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Table 63.7 Variation of anisotropic function with surface skin temperature 
for standard atmosphere; Atm Top = 30 km, spectral range 
- 5-50 km, A =3, RH * 75%, clear sky. 


Zenith 

Angle, 

e 

Surface Skin Temperature, K 


bbibo 

T 2 

T i=V 5 

T 2 

T z 

(C0 2 only) 



0 

15 

30 

45 

60 

70 

80 

1.074 

1.069 

1.053 

1.023 

0.974 

0.925 

0.847 

1.081 

1.075 

1.057 

1.025 

0.972 

0.917 

0.830 

1.079 

1.073 

1.056 

1.024 

0.972 

0.919 

0.835 

1.013 

1.012 

1.009 

1.004 

0.996 

0.988 

0.976 

1.083 

1.078 

1.059 

1.026 

0.971 

0.914 

0.823 

■ 

■ 

G 

0.227 

0.251 

0.244 

0.037 

0.260 

0.238 
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'V 


Table B3.8 


Variation of anisotropic function with 
concentration for wet and dry standard 
= 5-50 u, T. = 288.15 K, e. = 1.0, X = 


carbon dioxide 
atmosphere; Atm Top 
3, RH = 75X, cle^r sky. 


Zenith 

Angle, 

Carbon Dioxide Concentration 

Dry Atmosphere 

Standard Water Vapor 

e 

Zero 

Stand. 

Double 

Zero 

Stand. 

Double 

0 

15 

30 

45 

60 

70 

80 

1.013 

1.012 

1.009 

1.004 

0.995 

0.985 

0.967 

1.028 

1.026 

1.020 

1.009 

0.990 

0.970 

0.937 

1.027 
1.025 
1 = 019 
1.009 
0.990 
0.971 
0.938 

1.072 

1.067 

1.052 

1.023 

0.973 

0.918 

0.823 

1.079 

1.073 

1.056 

1.024 

0.972 

0.919 

0.835 

1.077 

1.072 

1.055 

1.024 

0.973 

0.921 

0.838 

G 

0.046 

0.091 

0.089 

0.249 

Q.244 

0.239 
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Table B3.9 Latitudinal variation of anisotropic function for 

climatological -average model atmospheres; Atm Top = 30 km, 
spectral range = 5-50 u, e = 1.0, A = 3, e r = 1.0, 

Z c = 6 km 


Zenith 

Angle, 

Climatological -Average Model Atmospheres 

Tropical 

Sub- Arc tic Winter 

Mid-Lat. Ave. 

e 

Cl ear 

Overcast 

Cl ear 

Overcast 

Cl ear 

Overcast 

0 

1.097 

1.062 

1.053 

1.025 

1.079 

1.040 

15 

1.090 

1.057 

1.049 

1.023 

1.073 

1.037 

30 

1.069 

1.044 

1.038 

1.018 

1.056 

1.029 

45 

1.030 

1.019 

1.017 

1.008 

1.024 

1.013 

60 

0.965- 

0.978 

0.981 

0.991 

0.972 

0.985 

70 

0.898 

0.933 

0.943 

0.971 

0.919 

0.957 

80 

0.789 

0.864 

0.882 

0.936 

0.835 

0.912 

G 

0.308 

— 

0.198 

0.171 

0.089 

0.244 

0.128 
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Table B3.10 Variation of anisotropic function for different model 
atmospheres; Atm Top = 30 km, spectral range 
= 5-50 u, e s s 1.0, X = 3, e c = 1.0, l Q - 6 km. 


Zenith 

Model Atmosphere 

- 

Angle, 

Sub -Arc. 

Summer, 

Mid-Lat. 

Winter, 

Mid-Lat. 

SuraneM 

e 

T $ = 288.45 K 

T $ = 272.59 K 

T $ = 296.22 K. 


Clear 

Overcast 

Cl ear 

Overcast 

Clear 

Overcast 

0 

1.071 

1.035 

1.065 

1.034 

1.083 

1.050 

15 

1.066 

1.033 

1.061 

1.032 

1.077 

1.047 

30 

1.050 

1.025 

1.047 

1.025 

1.059 

1.036 

45 

1.022 

1.011 

1.020 

1.011 

1.026 

1.016 

60 

0.975 

0.987 

0.977 

0.987 

0.971 

0.982 

70 

0.928 

0.963 

0.931 

0.961 

0.914 

0.947 

80 

0.855 

0.926 

0.859 

0.915 

0.824 

0.892 

G 

0.216 

0.109 

0.206 

0.119 

0.259 

0.158 
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Table B3.ll Latitudinal variation of anisotropic function for radiosonde- 
measured atmospheric models (Data 106 models, 9/29/58); Atm 
Top = 30 km, spectral range = 5-50 u, c $ - 1.0, clear sky. 


Zeni ch 
Angle, e 


Atmospheric Model 

- 

Havana, Cuba 

Nantucket, Mass. 

Thule, Greenland 

0 

1.0834 

1.0704 

1.0501 

15 

1.0776 

1,0655 

1.0467 

30 

1.0595 

1.0502 

1.0360 

45 

1.0261 

1.0220 

1.0159 

60 

0.9703 

0.9750 

0.9816 

70 

0.9125 

0.9266 

0.9448 

80 

0.8219 

0.8500 

0.8854 

G 

0.2615 

0.2204 

0.1647 
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Table B3.12 Seasonal variation of anisotropic function for radiosonde- 

measured atmospheric models (Data 106 model); Atm Top = 30 km, 
spectral range = 5-50 y, e s = 1.0, clear sky. 


— ■ ■ 

Zenith 

Tropical 

Mid-tat. 

Const. 

Sub-Arctic 

Angle, 

Keywest, 

Grand Jn., 

Denver, 

Eureka, 

e 

FI orida 

Colo. 

Colo. 

NWT 


08/01/58 

02/01/58 

08/01/58 

01/01/58 

08/01/58 

03/01/58 

0 

1.0835 

1.0815 

1.0718 

1.0565 

1.0561 

1.0322 

15 

1.0777 

1.0759 

1.0668 

1.0527 

1.0522 

1.0301 

30 

1.0594 

1.0581 

1.0512 

1.0406 

1.0398 

1.0234 

45 

1.0260 

1.0254 

1.0224 

1.0179 

1.0173 

1.0105 

60 

0.9706 

0.9711 

0.9746 

0.9794 

0.9806 

0.9875 

70 

0.9134 

0.9151 

0.925& 

0.9388 

0.9442 

0.9615 

80 

0.8226 

0.8278 

0.8489 

0.8749 

0.8905 

0.9150 

G 

0.2609 

0.2538 

0.2229 

0.1816 

0.1655 

0.1172 



APPENDIX B4 


ANISOTROPIC FUNCTIONS FOR SPECTRAL 
RANGE 5-200 y - Tables B4.1 - B4.16 


Table B4.1 Variation of anisotropic function with high-cloud cover 
fraction for standard atmosphere; spectral range 
= 5-200 u. Atm Top = 30 km, T = 288.15 * , e = 1.0, 

X = 3. e c = 1.0, Z r = 6 km, RH = 75%. 


Zenith 


Cloud Cover, % 

Angle, 




e 

0 

50 

100 

0 

1.0768 

1.0608 

1.0395 

15 

1.0714 

1.0566 

1.0368 

30 

1.0546 

1.0433 

1.0283 

45 

1.0238 

1.0189 

1.0125 

60 

0.9731 

0.9785 

0.9857 

70 

0.9212 

0.9367 

0.9575 

80 

0.8382 

0.8702 

0.9130 

G 

0.2385 

0.1906 

0.1265 


It 










Table B4.2 Variation of anisotropic function with cloud height for standard 
atmosphere; spectral range = 5-200 y, Atm Top = 30 km, 

T s = 228.15 K, e = 1.0, \ = 3, e = 1.0, KH = 75%. 


Zenith 


Cloud Height 

(Z c ), km 

- 


Angle, 







0 

50% Cloud 

lover 

100% Cloud 

Cover 


2 

6 

10 

2 

6 

10 

0 

1.0715 

1.0608 

1.0511 

1.0658 

1.0395 

1.0060 

15 

1.0665 

1.0566 

1.0475 

1.0612 

1.0368 

1.0056 

30 

1.0580 

1.0433 

1.0363 

1.0468 

1.0283 

1.0043 

45 

1.0221 

1.0189 

1.0158 

1.0204 

1.0125 

1.0019 

60 

0.9750 

0.9785 

0.9821 

0,9770 

0.9857 

0.9978 

70 

0.9268 

0.9367 

0.9475 

0.9329 

0.9575 

0.9936 

80 

0.8508 

0.8702 

0.8925 

0.8645 

0.9130 

0.9874 

G 

0.2207 

0.1906 

0.1586 

0.2013 

0.1265 

0.0186 
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Table B4.4 Variation of anisotropic function with surface relative 
hunldity for standard atmosphere; spectra 1 range = 5-2 
Atm Top = 30 km, T c = 288.15 K, = 1.0, X * 3, e„ = 


Zenith 


Angle, 

e 


Surface Relative Humidity (RH), % 


07 
07 
05 
02 
0.9737 
0.9229 
0.8422 


Cle-'r Sky 


.0768 

,.0714 

1.0546 

1.0238 

0.9731 

0.9212 

0.8382 


07 

07 

05 

02 


0.8725 

0.9194 

0.8346 


1.0384 

1.0358 

1.0276 

1.0122 

0.9859 

0.9579 

0.9123 


Overcast sky 


1.0395 

1.0368 

1.0276 

1.0125 

0.9857 

0.9575 

0.9130 


0.2331 0.2386 0.2438 0.1256 


0.1267 














Table B4.5 Variation of anisotropic function with water vapor scale-height 
for standard atmosphere; spectral range = 5-200 u. Atm Top 
= 30 km, T s = 288.15 K, e $ = 1.0, e c = 1.0, Z c = 6 km, RH = 75%. 


Zenith 

Angle, 

9 

Water Vapor Scale-Height Parameter, \ 

. 

Clear Sky 

Overcast Sky 


2 

3 

4 

2 

3 

4 

0 

1.0811 

1.0768 

1.0721 

■ 

1.0417 

HU 

1.0350 

15 

1.0754 

1.0714 

1.0671 

1 

1.0388 

msm 

1.0327 

30 

1.0575 

1.0546 

1.0514 

1 

1.0297 

1.0283 

1.0253 

45 

1.0250 


1.0225 


1.0129 

1.0125 

1.0113 

60 

0.9719 

0.9731 

0.9744 


0.9854 

C.9857 

0.9868 

70 

0.9177 

0.9212 

0.9247 


0.9578 

0.9575 

0.9601 

80 

0.8321 

0.8382 

0.8441 


0.9176 

0.9130 

0.9153 

G 

0.2490 

0.2385 

0.2280 

0.1241 

0.1265 

0.1197 

















Table B4.6 Variation of anisotropic function with surface emissivity for 
standard atmosphere; spectral range = 5-200 y, Atm Top 


= 30 km, T e = 288.15 K, A - 3, = e r = 1.0, 


z c = 


6 km, RH = 75%. 


Zenith 

Angle, 

e 


Surface Emissivity, 

e s 


Clear Sky 


50% Cloud 

Cover 


0.8 

0.9 

1.0 

0.8 

0.9 

1.0 

0 

1.0647 

1.0710 

1.0768 

1.0535 

1.0572 

1.0608 

15 

1.0602 

1.0660 

1.0714 

1.0498 

1.0532 

1.0566 

30 

1.0460 

1.0504 

1.0546 

1.0381 

1.0408 

1.0433 

45 

1.0200 

1.022Q 

1.0238 

1.0167 

1.C178 

1.0189 

60 

0.9775 

0.9762 

0.9731 

0.9811 

0.9798 

0.9785 

70 

0.9347 

0.9277 

0.9212 

0.9449 

0.9407 

0.9367 

80 

0.8681 

0.8527 

0.8382 

0.8881 

0.8790 

0.8702 

G 

0.1966 

0.2183 

0.2386 

0.1654 

0.1782 

0.1906 
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Table B4.7 Variation of anisotropic function high-cloud cover fraction 
for tropical atmosphere; spectral range * 5-200 y, Atm 
Top = 30 km, T = 302.59 K, e = 1.0, X = 3, r = 1.0, Z = 6 km, 

_ 7 not ^ ^ v v 


Zenith 


Angle, 

9 


12: 







Table B4.8 Variation of anisotropic function with cloud height for 
tropical atmosphere; spectral range = 5-200 u, Atm Top 
= 30. km, T * 302.59, e = 1.0, X = 3, e - 1.0, RH = 75%. 


Zenith 

Angle, 

0 

Cloud Height (Z ), km 
c 

50 Cloud Cover 

| 100% Cloud Cover j 


2 

6 

10 

2 

6 

10 

0 

1.0902 

1.0810 

1.0721 

1.0833 

1.0608 

1.0308 

15 

1.0840 

1.0754 

1.0671 

1.07761 

1.0567 

1.0287 

30 

1.0643 

1.0579 

1.0515 

1.0594 

1.0436 

1.0222 

45 

1.0282 

1.0254 

1.0227 

1.0260 

1.0192 

1.0099 

60 

0.9678 

0.9709 

0.9740 

0.9704 

0.9779 

0.9884 

70 

0.9047 

0.9136 

0.9227 

0.9124 

0.9342 

0.9650 

80 

0.8036 

0.8225 

0.8225 

0.8403 

0.8187 

0.8560 

G 

0.2866 

0.2585 

0.2318 

0.2646 

0.9158 

0.1050 
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Table B4.9 Variation of anisotropic function with high-cloud emissivity 
for tropical atmosphere; spectral range = 5-200 M, Atm Top 
= 30 km, T s = 302.59 K, e $ = 1.0, X = 3, RH = 75%. 


Zeni th 
Angle, 
e 

High Cloud (Z c =■ 10 km) 

Emissivity, e c 

50% Cloud Cover 

100% Cloud Cover 


0.5 

1.0 

0.5 

1.0 

0 

1.0858 

1.0721 

1.0721 

1.0308 

15 

1.0799 

1.0671 

1-0671 

1.0287 

30 

1.0612 

1.0515 

1.0515 

1.0222 

45 

1.0269 

1.0227 

1.0227 

1.0099 

60 

0.96921 

0.9740 

0.9740 

0.9884 

70 

0.9087 

0.9227 

0.9227 

0.9650 

80 

0.8120 

0.8403 

0.8403 

0.9258 

G 

0.2738 

0.2318 

0.2318 

0.1050 
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Table 84.10 Variation of anisotropic function with surface relative humidity 
for tropical atmosphere; special range = 5-200 y. Atm Top 
= 30 km, T s = 302.59 K, = 1.0, A = 3, RH = 75% 


Zenith 

Angle, 

e 

r 

Surface 

Relative Humidity (RH), % 



Clear Sky 

Overcast Sky 


50 

75 

100 

50 

75 

ICO 

0 

1.0917 

1.0967 

1.1009 

1.0588 

1.0608 

1.0621 

15 

1.0854 

1.0900 

1.0939 

1.0548 

1.0567 

1.0579 

30 

1.0654 

1.0690 

1.0720 

1.0423 

1.0436 

1.0445 

45 

1.0287 

1.0303 

1.0315 

1.0187 

1.0192 

1.0196 

60 

0.9671 

0.9654 

0.9641 

0.9784 

0.9779 

0.9776 

70 

0.9024 

0.8976 

0.8944 

0.9354 

0.9342 

0.9336 

80 

0.7972 

0.7894 

0.7859 

0.8660 

0.8650 

0.8645 

G 

0.2945 

0.3073 

0.3150 

0.1958 

0.1958 

u- 

0.1976 
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Table B4.ll Variation of anisotropic function with water vapor scale height 
parameter for tropicaL atmosphere; spectral range = 5-200 u, 

Atm Top = 30 km, T s = 302.59 K, e $ = 1.0, Z c = 6 km, RH = 75% 


Zenith 

Angle, 

0 

Water Vapor Scale-He 

ight Parameter, A 


Clear Sky 


Overcast 

>ky 


2 

3 

4 

2 

3 

4 

0 

1.1103 


1.0861 

1.0671 

1.0608 

1.0522 

15 

1.1026 

1.0900 

1.0802 

1.0624 

1.0567 

1.0488 

30 

1.0784 

1.0690 

1.0615 

1.0476 

1.0436 

1.0377 

45 

1.0342 

1.0303 

1.0271 

1.0207 

1.0192 

1.0167 

60 

0.9614 

0.9654 

0.9689 

0.9768 

0.9779 

0.9805 

70 

0.8878 

0.8976 

0.9074 

0.9332 

0.9342 

0.9410 

80 

0.7782 

0.7894 

0.8066 

0.8702 

0.8650 

0.8747 

G 

0.3321 

0.3073 

0.2795 

0.1969 

0.1958 

0.1775 
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Table B4.12 Variation of anisotropic function with surface emissivity 
for tropical atmosphere; spectral range = 5-200 y, Atm Top 
= 30 km, T s = 302.59 K, A = 3, e= 1.0, l = 6 km, Rli = 75%. 


Zenith 

Angle, 

9 


Surface Emissivity, e c 

- 


Clear Sky 

50% Cloud 

Cover 


0.8 

0.9 

1.0 

0.8 

0.9 

1.0 

0 

1.0787 

1.0879 

1.3967 

1.0707 

1.0759 

1.0810 

15 

1.0732 

1.0818 

1.0900 

0.0658 

1.0707 

1.0754 

30 

1.0561 

1.0627 

1.0690 

1.0505 

1.0542 

1.0579 

45 

1.0246 

1.0275 

1.0303 

1.0222 

1.0238 

1.0254 

60 

0.9717 

0.9685 

0.9654 

0.9745 

0.9727 

0.9709 

70 

0.9154 

0.9063 

0.8976 

0.9238 

0.9187 

0.9136 

80 

0.8197 

0.8042 

0.7894 

0.8400 

0.8311 

0.8225 

G 

0,2590 

0.2837 

0.3073 

0.2307 

0.2448 

0.2585 


126 




Table B4.13. Variation of anisotropic function with cloud cover fraction 

for subarctic-winter atmosphere; spectral range = 5-200 y, Atm 
Top = 30 km, T s = 257.28 K, e s = 1.0, X = 3, e c = 1.0, Z c = 6 
km, RH = 75%. 


Zenith 

Angle, 

e 

Cloud Cover, % 

0 

50 

100 

0 

1.0493 

1.0425 

1.0225 

15 

1.0459 

1.0396 


30 

i.0354 

1.0270 

1.0163 

45 

1.0156 

1.0120 

1.0073 

60 

0.9819 

0.9861 

0.9914 

70 

0.9461 

0.9582 

0.9736 

80 

0.8871 

0.9115 

0.9423 

G 

0.1622 

0.1310 

0.0802 














Table 84.14 Variation of anisotropic function with cloud height for sub- 
arctic winter atmosphere; spectral range = 5-200 y, Atm Top 
= 30 km, T s = 257.28 K, e $ = 1.0, \ = 3, e c 1.0, RH = 75%. 


Zenith 

Angle, 

e 

Cloud Height 

;z c ), km 

- 


50% Cloud Cc 

)ver 

100% Cloud Cover 


2 

6 

10 

2 

6 

■B 

0 

1.0485 

1.0425 

Bff * S t 

1.0477 

1.0225 

wmm 

15 

1.0452 

1.0396 

Bln 5: 9 

1.0445 

1.0210 

Hni! 9 

30 

1.0348 

1.G270 

mm c \ 9 

1.0343 

1.0163 

BSI 9 

45 

1.0154 

1.0120 

i . 

1.0152 

1.0073 

BiS» 

60 

0.9822 

0.9861 

§11 IK : : 9 

0.9825 

0.9914 

1 ? 9 

70 

0.9468 

0.9582 

0.9661 

0.9475 

0.9736 

■9?*' 9 

80 

0.8884 

0.9115 

0.9288 

0.8898 

0.9423 

0.9934 

G 

0.1601 

0.1310 

0.1021 

— 

0.1598 

0.0802 

0.0089 
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Table B4.15 Variation of anisotropic function with high-cloud emissivity 
for sub-arctic winter atmosphere; spectral range * 5-200 y, 
Abu Top = 30 km, T = 257.28 K, e - 1.0, X = 3, RH = 75%. 


Zeni th 
Angle, 
9 

High Cloud (Z c = 10 km) 

Emissivity, e c 

50% Cloud Cover 

100% Cloud Cover 


0.5 

1.0 

0.5 

1.0 

0 

1.0410 

1.0309 

1.0309 

1.0023 

15 

1.0383 

1.0288 

1.0288 

1.0022 

30 

1.0294 

1.0222 

1.0222 

1.0017 

40 

1.0130 

1.0098 

1.0098 

1.0008 

60 

0.9850 

0.9887 

0.9887 

0.9991 

70 

0.9551 

0.9661 

0.9661 

0.9971 

80 

0.9050 

0.9288 

0.9288 

0.9934 

G 

0.1351 

0.1021 

0.1021 

0.0089 
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Table B4.16 Variation of anisotropic function with surface emissivity for 
sub-arctic winter atmosphere; spectral range ~ 5-200 y. Atm 
Top = 30 km, T c = 257.28 K, \ = 3, e r 1.0 Z r = 1.0 L =6 km, 
RH = 75%. 5 


Surface Emissivity, e, 


Clear Sky 


50% Cloud Cover 


0.8 0.9 


1.0261 

1.0118 

0.9860 

0.9569 

0.9066 


0.8965 


1.0354 1.0216 


RI 


0.9873 

0.9614 

0.9172 


0.1293 0.1463 


0.1622 0.106610.1165 
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APPENDIX Cl 

SYMBOLS USED IN THE COMPUTER PROGRAM "FILAUP 


APPENDIX Cl 

SYMBOLS USED IN THE COMPUTER PROGRAM "FILAUPG 


AISFUN 

Anisotropic function 

ALA 

Altitude dependent average width of the lines of 
a molecule. 

'•vLB 

Altitude dependent individual width of the lines of 
a molecule. 

AVSI 

Average value of intensity 'or the lines in one 
decade in an interval. 

AVS1 


AVS2 


AVS3 

Average intensities of N 2 0, CH 4 , C0 2 , H 2 0 and 0 3 

AVS4 


AVS5 


CC 

Cloud cover 

CEMI 

Cloud emissivity 

DEL 

width of an interval 

EMG 

Surface emittance 

FLUX 

'Jpwelling flux at „he top of the atmosphere, W/M 2 

FRL 

Lower frequency limit of the range. 

IC 

Integer which is equal to zero if effect of cloud 
is not considered and one if it is considered. 

11, 12 

Integer for each of the five gases considered is 

13, 14 

equal to 1. It is equal to zero if the effect of a 

15 

is neglected. 

JD 

Number of adjacent intervals on both sides of an 
interval from which contribution is taken into 
account. 

KR 

Number of intervals in the band 

LCB 

Cloud height 

LT 

Height of the atmosphere 
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NL 

Number of layers 

NS I 

NS1, NS 2 
NS3, NS4 
NS5 

Number of lines in a decade within an interval 
Number of lines in an interval. 

PL 

Optical path length at the frequency under 
consideration. 

PLG 

Planck function of atmosphere evaluated at the 
average temperature of the layer. 

PNTP 

Pressure at NTP, bar. 

PRG 

Average pressure of a layer. 

PSC 

Planck function of the cloud evalulated at the 
average temperature of the cloud 

PSG 

Planck function of the surface evaluated at the 
surface temperature. 

RADC1 , RADC2 
RADC3.RADC4 
RA0C5.RADC6 

Radiances in each of the six spectral ranges 
considered. Watts M-2 Sr-1. 

RDAG 

Radiance emitted by the atmosphere, watts M-2 SR-1 

R DIG 

Total Radiance at the top of the atmosphere, watts 
M-2 Sr-1. 

RDSC 

Radiance emitted by the cloud. Watts M-2 Sr-1. 

RDSG 

Radiance emitted by the surface, watts M-2 Sr-1. 

TAG 

Combined transmittance of all interfering gases. 

TEG 

Average temperature of a layer, K. 

TEMPG 

Surface temperature, K. 

TEMR 

Reference temperature for the line parameters, K. 

TH 

Thickness of a layer the atmospheric layer, km. 

TNTP 

Temperature of NTP, K. 

TRAC 

Cloud transmittance 

VLIMDF 

Limb darkening function 


133 


NL 

Number of layers 

NSI 

NSI, NS 2 
NS3, NS4 
NS5 

Number of lines in a decade within an interval 
Number of lines in an interval. 

PL 

Optical path length at the frequency under 
consideration. 

PLG 

Planck function of atmosphere evaluated at the 
average temperature of the layer. 

PNTP 

Pressure at NTP, bar. 

PRG 

Average pressure of a layer. 

PSC 

Planck function of the cloud evalulated at the 
average temperature of the cloud 

PSG 

Planck function of the surface evaluated at the 
surface temperature. 

RADC1 , RADC2 
RADC3,RADC4 
RADC5.RADC6 

Radiances in each of the six spectral ranges 
considered. Watts M-2 Sr-1. 

RDAG 

Radiance emitted by the atmosphere, watts M-2 SR-1 

RDIG 

Total Radiance at the top of the atmosphere, watts 
M-2 Sr-1. 

RDSC 

Radiance emitted by the cloud. Watts M-2 Sr-1. 

RDSG 

Radiance emitted by the surface, watts M-2 Sr-1. 

TAG 

Combined transmittance of all interfering gases. 

TEG 

Average temperature of a layer, K. 

TEMPG 

Surface temperature, K. 

TEMR 

Reference temperature for the line parameters, K. 

TH 

Thickness of a layer the atmospheric layer, krn. 

TNTP 

Temperature of NTP, K. 

TRAC 

Cloud transmittance 

VLIMDF 

Limb darkening function 
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APPENDIX C2 

LISTINGS OF COMPUTER PROGRAM f, FILAUPG 


o o ooooooo 


APPENDIX C2 



LISTINGS OF COMPUTER PROGRAM "FILAUFG" 


PROGRAM FILAUPG ( INPUT , OUTPUT , TAPE5 - INPUT , TAPE6= OUTPUT , TAPE2 ) 

C 

C 

C THIS RADIATIVE TRANSFER PROGRAM USES QUASI-RANDOM BAND MODEL. 

C THE LONGWAVE REGION (50 - 2000CM- 1 ) IS DIVIDED INTO 195 INTERVALS 
C OF 1 0CM__ 1 EACH. THE ATMOSPHERE T30KM) IS DIVIDED INTO 15 LAYERS OF 
C THICKNESSES 1 ,1,1,1 ,1 ,1,1 ,1 f 1 ,1 ,1 , 4 , 5 , 5 , 5KM . PURE HTATIONAL AND 
C 6 . 3M ( MICRON ) BANDS AND CONTINUUM BAND OF WATER VAPOR IN 8-28. 5M 
C REGION, 15M CARBON DIOXIDE BAND, 9 ♦ 6M OZONE BAND, 7-8M NITRUS OXIDE 

C BAND AND 7.7M METHANE BAND ARE CONSIDERED. SINCE NUMEER OF LINES 

IN SOME BANDS IS VERY LARGE, PRECOMPUTED VALUES OF AVERAGED 
INTENSITY AND NUMBER OF LINES FOR EVERY INTERVAL AND INTENSITY 
DECADE ARE ISED IN THIS PROGRAM.’ TRANSMITTANCE COMPUTATION FOR 
EACH BAND IS CARRIEDOUT IN SUBROUTINE TRANS WHICH USES A SIMPLE 
INTEGRATION SCHEME EXPLAINED IN REPORT 75-T14. 


INTEGER X , TH 

DIMENSION PLG(1 95 , 1 5 ) ,RDSG(7) , RDAG ( 7 ) , RDSC (7 ) 

DIMENSION VLIMDF (7,4) 

DIMENSION SUG(5 ) ,SPUG(5 ) ,STUG(5) ,ULG(l 5 ,5) ,PULG(15 , 5 ) , TULG ( 1 5 , 5 ) 
DIMENSION PR(30) ,TE(30) ,VMR(30,5) , UG ( 1 5 , 5 ) ,PRG(1 5,5 ) , TEG(15 ,5 ) 
DIMENSION PREG ( 1 5,5), TEMG ( 1 5 , 5 ) , ALA (5 ) , ALG ( 1 5 , 5 ) , TH ( 1 5 ) , ALT ( 30 ) 
DIMENSION AVS1 ( 5 , 24 ) , AVS2 ( 5 , 2 9 ) , AVS3 ( 5 , 25 ) , AVS4 ( 5 , 1 95 ) , AVS5 (5 , 23 ) 
DIMENSION NS1 (5 ,24) ,NS2(5,29) ,NS3(5 ,25) ,NS4(5 , 195) ,NS5(5,23) 
DIMENSION TRG(1 95,16) ,PPG( 1 6 ) ,FAG( 1 6 ) 

DIMENSION AG (90) ,ZX(7 ) , ZEN (7) , RADG ( 7 ) ,TAG(7 , 195,16) 

DIMENSION RDI G ( 7 ) ,RDG(1 95) ,FXG(195) 

DIMENSION PNUM ( 1 95 ) ,EEX( 1 95) ,PSC(l 95 ) ,PSC(l 95 ) 

DIMENSION A IS FUN (7 ,5 ,2) , RADCC (7,5,2) , F LUX ( 5 , 2 ) , CC ( 5 ) 

DIMENSION TR AC ( 1 1 ) ,CEMl(l 1 ) ,RADNCL(l95) 

DIMENSION RADC1 (7,14) ,RADC2(7 , 1 4) ,RADC3(7 , 14) ,RADC4(7 , 1 4 ) 
DIMENSION RADC5(7,14) ,RADC6(7,14) ,RADZ(7,195) 

COMMON /TRANE/ AVSI (5 , 1 95) , NSI ( 5 , 1 95 ) , FRC ( 1 95 ) ,TRA(195) , 

1X1 (26) ,T1 (26) ,X2(21 ) ,T2(21 ) , DELA , JD , PI 

ZENITH ANGLE 
C 

DATA ZEN/O. , 1 5 • , 30 . ,45 • , 60 . ,70 . ,80 . / 

C 
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ORIGINAL p.A'%£ rr* 

OF POOR QUA* ITY 

C CLOUD COVER FRACTION T 

C 

DATA CC/O. 0,0. 3, 0.5, 0.7,1 .0/ 

DATA PNTP,TNTP, TEHR/l .00,273-1 5,296./ 

DATA ALA/O. 079, 0.05 4, 0.066, 0.063, 0.1 10/ 

DATA UG/75*0 . / 

DATA FRL1 ,FRL2 ,FRL3 ,FRL4 ,FRL5/1 1 1 0 . , 1 220 . , 5 50 . , 50 . ,950./ 

DATA DEL, DELA,JD, PI/10. ,5-, 4, 5- 1459/ 

DATA KR1 ,KR2,KR3,KR4,KR5 , KRT/24 , 29 , 25 , 1 95 , 23 , 1 95/ 

DATA KB1 ,02,03 , KB4 , KB5 /I 07 , 1 18,51 .1 ,91/ 

DATA KE1 ,KE2,KE3,KE4,KE5/150,146,75,195,115/ 

DATA FRLC,KC,KBC,KEC/350. ,90,31 ,120/ 

DATA A, B, BETA/1 . 25E-22 , 2 . 34 E-1 9 , 8 . 30E-O3/ 

REWIND 2 
C 

C VOLUME MIXING RATIOS N20,CH4,C02 

C 

DO 4 L-1 ,30 
VMR ( L , 1 )-0 .28 
VMR (L , 2 ) * 1 .6 
VMR (L , 3 ) “330 . 

4 CONTINUE 

DATA IC ,RE/1 ,20 ./ 

C 

C ATMOSPHERE IS DIVIDED INTO 15 LAYERS 

C 

DATA TH/1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,4, 5, 5,5/ 

C 

C PRESSURE AND TEMPERATURE PROFILE. VOLUME MIXING RATIOS 
C OF WATER VAPOR AND OZONE 
C 

READ(5 ,*) (PR(L) ,TE(L) , VMR ( L , 4 ) , VMR ( L , 5 ) , L- 1 , 30 ) 

C 

C FIVE GASES CONSIDERED ARE N20 , CK4 , C02 , H20 , 03 

C 

DATA II ,12,13.14,15/1 ,1 ,1 ,1 , 1/ 

C 

C NUMBER OF LAYERS CONSIDERED AND THE HEIGHT OF THE ATMOSPHERE 

C 

NL= 1 0 IF ATMOSPHERE TOP IS AT 1 OKM 
NL* 1 3 IF ATMOSPHERE TOP IS AT 20KM 

DATA NL , LT/1 5 ,30/ 

DATA XI /O .0 ,0 .001 ,0.001 5 ,0.002,0.003,0.004,0.005 ,0.006 , 

10.008,0.01 ,0.015 ,0-02,0.03,0.04,0.05,0.06,0.08,0.10, 

20.15.0. 20.0.30.0.40.0.50.0.60.0.80.1.00/ 

DATA NL , LT/1 5 , 30/ 

DATA T1 /O .0006 ,0.0006 ,0.0006,0.0007,0.001 ,0.001,0.001, 

10.0015 .0. 002.0.003 .0.005 .0.008.0.01 . 0. 01, C. 01, 0.015, 

20.02.0. 03.0.05.0.08.0.1 ,0.1 ,0.1 ,0.15,0.2,0.1/ 

DATA X2/-1 .0,-0. 9, -0.8, -0.7, -0.6 ,~0.5 ,-0.4 ,-0.3 ,-0.2, 

1-0.1 ,0.0, 0.1 ,0.2, 0.3, 0.4, 0-5, 0.6, 0.7, 0.8, 0.9,1 .0/ 

DATA T2/0. 1,0. 4, 0.2, 0.4, 0.2, 0.4, 0.2, 0.4, 0.2, 

10.4.0. 2.0.4.0.2.0.4.0.2.0.4.0.2.0.4.0.2.0.4.0.1/ 

C 

C AVERAGE LINE INTENSITIES AND NUMBER OF LINES IN EACH INTERVAL 
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C 


100 

C 

C 

C 

C 


READ(2 ,*) ( ( AVS1 ( I , K ) , I- 1 f 5 ) , K» 1 , KR1 ) 
READ(2 , *) ( (NS1 (l,K) ,1-1 ,5) ,K-1 t KR1 ) 
READ (2,*) ((AVS2(I,K) , 1-1 ,5) ,K-1 ,KR2) 
REA T> ( 2 , * ) ((NS2(I,K) ,1*1 ,5) ,K»1 ,KR2) 
READ(2 ,*) ((AVS3(I,K) ,1*1 ,5),K-1 ,KR3) 
READ (2,*) ((NS3(I,K) ,1-1 ,5) ,K*1 ,KR3) 
RE AD (2,*) ( ( AVS4 ( I , K ) , I* 1 ,5),K-1 ,KR4) 
READ (2,*) ((NS4(I,K),I-1 ,5) ,K-1 ,KR4) 
READ(2 ,*) ((AVS5(I,K) ,1-1 ,5),K-1 ,KR5) 
READ(2 ,*) ((HS5(I,K) ,1-1 ,5) ,K-1 ,KR5) 
DO 100 L* 1 ,30 
ALT (L ) -L-0.5 
C-0.1*TNTP/PNTP 


ORIGINAL PASL' & 
OF POOR QUALITY 


AVERAGE TEMPERATURE AND PRESSURE FOR ALL THE LAYERS OF 
OF ATMOSPHERE ARE COMPUTED USING CURTIS GODSON APPROXIMATION 


DO 101 N-1 ,5 
L2-LT 

SUG(N)-SPUG(N)=STUG(N)=0. 
DO 102 M-1 ,NL 
LA-NL+1 -M 


ulg(la,n)=pulg(la,n)-tulg(la,n)-o. 

LI -L2-TH( LA) +1 

DO 103 L-L1 ,L2 
DU=C*PR(L)*VMR(L,N)/TE(L) 

ULG(LA,N)=ULG(LA, N)+DU 
PULG(LA,N)=PULG(LA,N)+PR(L)*DU 
TULGCLA,N)-TULG(LA,N)+TE(L)*DU 
SUG (n)-SUG(n)+DU 
spug(n)-spug(n)+pr(l)*du 
103 stug(n)-stug(n)+te(l)*du 
ug(la,n)=sug(n) 

prg(la,n)=pulg(la,n)/ulg(la,n) 
teg(la,n)-tulg(la,n)/ulg(la,n) 
preg(la,n)=spug(n)/sug(n) 
temg(la, n)-stug(n)/sug(n) 

alg(la,n)=ala(n)*(sqrt(temr/temg(la,n) ) )*PREG(LA,N)/PNTP 
102 L2-L1-1 

101 CONTINUE 

c 

TRANSMITTANCE FOR THE N2C BAND IS COMPUTED FOR ALL THE 15 
LAYERS. TRG IS THE FINAL TRANSMITTANCE, SO INDIVIDUAL 
TRANSMITTANCES ARE MULTIPLIED INTO THIS ARRAY 
C 

LD-NL+1 

DO 110 LA-1 , LD 
DO 110 K-1 , KRT 
110 TRG ( K , LA ) -1 . 

IF ( 1 1 . LT . 1 ) GO TO 163 
DO 111 K-1 ,KR1 
DO 111 1-1,5 
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1 1 1 


113 
1 1 2 
C 
C 
C 

163 


116 


118 

117 

C 

C 

C 

164 


FRC(K)-FRL1 +(2*K-1 ) *DELA 
AVSI ( I ^ -AVS 1 ( I , K ) 

NSI ( I , K , =NS1 ( I , K ) 

CONTINUE 

DO 1 1 2 LA-1 ,NL 

CALL TRANS (ALG( LA, 1 ) , UG ( LA , 1 ) , KR1 ) 
DO 1 1 3 K-KB1 , KE1 

TRG(K,LA)-TRG(K,LA)*TRA(K-KB1 +1 ) 
CONTINUE 


OF POOR Q UAiLr.ir 


TRANSMITTANCE FOR CH4 BAND 

IF ( 12 . LT . 1 ) GO T0164 
DO 1 1 6 K= 1 , KR2 
FRC(K)-FRL2+(2*K-1 )*DELA 
DO 116 1=1 ,5 
AVSI(I,K)=AVS2(I,K) 

NSI(I ,K)=NS2(I,K) 

CONTINUE 

DO 1 17 LA* 1 , NL 

CALL TRANS(ALG(LA,2) ,UG(LA,2) ,KR2) 
DO 118 K-KB2 , KE2 

TRG(K,LA)=TRG(K,LA)*TRA(K-KB2+1 ) 
CONTINUE 

TRANSMITTANCE FOR C02 BAND 

I F ( 13 • LT . 1 ) GO TO 165 
DO 121 K= 1 ,KR3 
FRC(K)=FRL3+(2*K-1 )*DELA 


DO 121 1 = 1 ,5 
AVSI(I ,K)-AVS3(I ,K) 

NSI(I ,K)=NS3(I,K) 

121 CONTINUE 

DO 1 22 L A = 1 , NL 

CALL TRANS ( ALG (LA , 3 ) , UG ( LA , 3 ) , KR3 ) 

DO 123 K=KB3,KE3 

123 TRG(K,LA)=TRG(K,I<A)*TRA(K-KB3 + 1 ) 

122 CONTtNUE 
C 

C TRANSMITTANCE FOR WATER VAPOR FOR 6.3M AND PURE 

C ROTATIONAL BAND 

C 

165 IF (I4.LT.1 ) GO TO 166 
DO 1 26 K* 1 ,KR4 
FRC(K)=FRL4+(2*K-1 ) *DELA 
DO 126 1=1,5 
AVSI(I,K)=AVS4(I,K) 

NSI(I,K)-NS4(I,K) 

126 CONTINUE 

DO 1 27 L A= 1 , NL 

CALL TRANS (ALG (LA, 4) ,UG(LA,4) ,KR4) 

DO 128 K=KB4,KE4 

1 28 TRG(K,LA)=TRG(K,LA)*TRA(K-KB4 + 1 ) 

127 CONTINUE 139 


ORIGiftfciL. [.>j 

OF POOR QUALITY 

C 

C TRANSMITTANCE IS COMPTED CORRESPONDING TO WATER VAPOR 

C CONTINUUM ABSORPTION INTHE 350-1250CM 1 

C ” 

DO 131 K-1 , KC 

131 FRC (K) -FRLC+ (2*K- 1 )*DELA 
DO 1 32 LA-1 , NL 

PPG(LA)=PREG(LA,4)*UG(LA,4)*VMR(l ,3)*1 • E-06/UG ( LA , 3 ) 

132 FAG(LA)=(PPG(LA)*EXP(6 .08*( ( TEMR/ TEMG ( LA , 4 ) ) 

1 -1 . ) ) + 0.002*(PREG(LA,4)-PPG(LA)) ) 

DO 134 LA-1 , NL 
DO 135 K-1 , KC 

135 AG(K)=2.69E+1 9*(A+B*EXP(-BETA*FRC(K) ) )*FAG(LA) 

DO 137 K-KBC , KEC 

137 TRG(K,LA)-TRG(K, L A )*EXP ( -UG ( LA , 4 ) *AG (K-KBC+ 1 )) 

C 

C TRANSMITTANCE FOR 03 BAND 

C 

134 CONTINUE 

1 06 IF (I5.LT.1) GO TO 167 
DO 142 K-1 ,KR5 
FRC(K)=FRL5+(2*K-1 )*DELA 
DO 142 1=1 ,5 
AVSI(I ,K)=AVS5(I ,K) 

NSI(I,K)-NS5(I,K) 

142 CONTINUE 

DO 145 LA-1 , NL 

CALL TRANS (ALG( LA , 5) ,UG(LA,5) ,KR5 ) 

DO 144 K-KB5 , KE5 

1 44 TRG(K,LA)=TRG(K,LA)*TRA(K-KB5+1 ) 

143 CONTINUE 
167 CONTINUE 

CONS-1 8. *6 .625E-03 
CNST-6. 625*0. 3/1 -38 
DO 147 M-1 ,7 

ZX ( M ) - 1 . /COS ( ZEN (M)/57. 29578) 

147 CONTINUE 


C 

C 

C PLANCK FUNCTIONS FOR EACH OF 195 INTERVALS AND THE 15 

C LAYERS OF ATMOSPERE ARE COMPUTED AND STORED 

C 
C 

DO 31 1 K-1 , KRT 
FRC(K)«50.+(2*K-1 ) *DEL A 

PNUM(K)=DEL*C0NS*FRC(K)*FRC(K)*FRC(K)*1 . E-07 
EEX(K)-CNST*FRC(K) 

DO 149 M-1 ,7 
DO 1 50 LA-1 , NL 
C 
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TRANSMITTANCES IN SOME OF THE INTERVALS ARE VERY LOW 
BECAUSE OF STRONC, ABSORPTION. TO AVIOD UNDERFLOW ERROR 
THEY ARE EQUATED TO ZERO. 

IF ( TRG (K , LA) . LT . 1 .E-50) GO TO 200 
TAG(M,K,LA)=TRG(K,1A)**ZX(M) 

GO TO 150 

200 TAG (M,K , LA) =0 . 

150 CONTINUE 

TAG ( M , K , NL+ 1 )«1 . 

149 CONTINUE 

DO 152 LA-1 ,5 

152 PLG(K,LA)*PNUM(K)/(EXP(EEX(K)/TEG(LA t 4))-1 . ) 

IF(NL.LE.5) GO TO 511 

DO 153 LA-6, NL 

153 PLG(K,LA)-PNUM(K)/(EXP(EEX(K)/TEG(LA,3) )-1 .) 

311 CONTINUE 

SURFACE TEMPERATURE 

TEMPG-302 . 59 
DO 213 KTE-1 ,3 
TEMPG-TEMPG+5 • 

DO 315 K-1 , KRT 

313 psg(k)=pnum(k)/(exp(eex(k)/tempg)-i . ) 

HEIGHT OF THE CLOUD BASE CAN BE CHANGED BY CHANGING LCB HERE 
LCB= -2 

DO 212 NK-1 ,3 
LCB-LCB+4 
LCT-LC B+ 1 
C 

C AVERAE CLOUD TEMPERATURE 

C 

TEMPC-0.5*(TE(LCB)+TE(LCB+1 ) ) 

DO 312 K-1 ,KRT 

312 PSC(K)-PNUM(K)/(EXP(EEX(K) /TEMPO )-1 . ) 

c 

C SURFACE EMISSIVITY 

C 

EMG-0.70 
DO 21 1 NE-1 ,3 
EMG-EMG+0.1 

209 F0RMAT("1 SURFACE EMISSIVITY - " , F6 . 2 , // , " 

1 SURFACE TEMPERATURE - M F6.2,//," 

2HEIGHT OF THE ATMOSPHERE - " , 12 , "KM" , // , " 

3CL0UD TOP HEIGHT-", 12, //,” 

4CL0UD EMISSIVITY =",F5.3,//) 
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DO 214 M- 1 ,7 

radg(m)*o. 

DO 214 X-1 ,6 

RADC1 (M,X)-RADC2(M,X)-RADC3(M,X)-RADC4(M,X)«RADC5(M,X) 

1 =RADC6(M,X)«0. 

214 CONTINUE 
FLG-O. 

CLEAR ATMOSPHERE UPWELLING RADINCES AND FLUXES FOR 
ALL ZEITH ANGLES ARE COMPUTED 

DO 148 K-1 , KRT 
DO 154 M-1 ,7 

UPWELLING RADIANCE DUE TO SURFACE EMISSION 
RDSG(M)-EMG*PSG(K)*TAG(M,K, 1 ) 

rdag(m) *0 . 

DO 155 LA-1 , NL 

UPWELLING RADIANCE DUE TO ATMOSPHERIC EMISSION 
1 55 RDAG{M)=RDAG(M)+PLG(K,LA)*(TAG(M,K,LA+1 )-TAG(M,K, LA)) 

rdig(m)-rdsg(m)+rdag(m) 

radg(m)-radg(m)+rdig(m) 

154 CONTINUE 

COMPUTES UPWELLING FLUX 3Y EXPONENTIATING TRANSMITTANCE 
BY THE DIFFUSIVITY FACTOR 1.66 AND MULTIPLYING PLANCK 
FUNCTION TIMES PI 

FLSG-EMG*PI*PSG(K)*(TRG(K, 1 )**1 .66) 

FLAG-0. 

DO 1 57 LA-1 , NL 

1 57 FLAG-FLAG+PI*PLG(K,LA)*( (TRG(K,LA+1 )**1 .66) 

1 -(TRG(K,LA)**1 .66) ) 

FLXG-FLSG+FLAG 
FLG-FLG+FLXG 
RDG (K )-RDIG ( 1 ) 

FXG(K)-FLXG 
148 CONTINUE 

DO 1 59 X-1 ,2 

ONLY TWO DIFFERENT CLOUD EMISSIVTIES ARE USED IN THIS 
STUDY. DIMENSION OF X SHOULD BE CHANGED IF IT IS REQUIRED 
TO VARY MORE THAN TWO 

FLC-O. 

CLOUD EMISSI VITY 

ONLY THE HIGH CLOUD HAS TWO DIFFERNT VALUES 
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IF(LCB.LE.6.AND.X.EQ.2) GO TO 159 
IF(LCB.LE.6.AND.X.EQ. 1 ) CEMI (X ) ■ 1 .0 
IF(LCB.EQ.10.AND.X.EQ.1 ) CEMl(X)-0.5 
IF(LCB.EQ. 10c AND . X . EQ . 2 ) CEMI(X)-1 .0 


CLOUD TRANSMITTANCE IS EVALUATED 

TRAC (X ) - 1 .-CEMI(X) 

DO 210 M- 1 ,7 
DO 210 K=* 1 , KRT 
210 RADZ ( M , K ) *0 . 

IF(IC.EQ.O) GO TO 601 
DO 51 5 K- 1 , KRT 

UPWELLING FLUX FROM CLOUD TOP 

FLSC-CEMI(X)*PI*PSC(K)*(TRG(K,LCT)**1 .66 ) 
FLSG=EMG*PI*PSG(K)*(TRG(K, 1 )**1 .66) 

FLAC-O. 

UPWELLING FLUX FROM THE ATMSPHERIC LAYERS ABOVE THE CLOUD 
DO 222 LA-LCT , NL 

222 FLAC*FLAC+PI*PLG(K,LA)*((TRG(K,LA+1 )**1 .66) 

1 -(TRG(K f LA'**1 .66) ) 

FLXC-FLSC+FLAC 

FLC-FLC+FLXC 

IF(TRAC(X) .EQ.O.O) GO TO 3909 
FLAT-0. 

IF THE CLOUD EMISSIVITY IS LESS THAN 1 f UPWELLING FLUX 
COMPONENT BELOW THE CLOUD LAYER IS EVALUATED 

DO 223 LA=1 ,LCB 

2 23 FLAT«FLAT+PI*PLG(K,LA)*( (TRG ( K , LA+ 1 ) **1 . 66 ) 

1 -(TRG(K, LA)**1 .66) ) 

FLAB- ( FLAT +FLSG ) *TRAC (X ) 

FLC-FLC + FL/.B 
909 DO 158 M-1 ,7 

UPWELLING RADIANCE FROM CLOUD TOP 

RDSC(M)-CEMI(X)*PSC(K)*TAG(M,K, LCT) 

UPWELLING RADIANCE FROM THE SURFACE AND THE ATMOSPHERE 
BELOW THE CLOUD 
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RDSG(M)-EMG*PSG(K)*TAG(M,K, 1 ) 

RADNCL (K ) -0 . 

DO 409 LA-1 , LCB 

TERM2*PLG(K,LA) # (TAG(M,K,LA+1 )-tag(m,k,la")) 

409 RADNCL(K)-RADNCL(K)+TERM2 

RADNCL (K) -TRAC (X )* (RDSG (M) +RADNCL (K) ) 

C EMISSION PROM CLOUD 

RADNCL (K) -RADNCL (K)+RDSC(M) 

U.PWELLING RADIANCE FROM THE ATMOSPHERE ABOVE THE CLOUD 
DO 408 LA-LCT , NL 

term4*plg(k,la) # (tag(m,k,la+i)-tag(m,k,la)) 

408 RADNCL(K)-RADNCL(K)+TERM4 


UPVELLING RADIANCE AT THE TOP OF THE ATMOSPHERE FOR A 
PARTICULAR ZENITH ANGLE AND A SUB-INTERVAL 

RADZ(M,K)-RADNCL(K) 

58 CONTINUE 
15 CONTINUE 
01 CONTINUE 


DO 190 M-1 ,7 

TOTAL UPVELLING RADIANCE IN THE SPECTRAL REGION 5-1 OM 
DO 180 K-96,195 

RADC1 ( M , X) -RADC1 ( M , X ) +R ADZ (M , K ) 

80 CONTINUE 

TOTAL UPVELLING RADIANCE IN THE SPECTRAL REGION 10-20M 
DO 181 K-46,95 

RADC2(M,X)-RADC2(M,X)+RADZ(M,K) 

81 CONTINUE 

TOTAL UPVELLING RADIANCE IN THE SPECTRAL REGION 5-20M 
RADC3 (M,X)-RADC1 ( M , X) +RADC2 ( M , X ) 

TOTAL UPVELLING RADIANCE IN THE VINDOV REGION 10.5-12-5 
DO 182 K-76,90 

RADC4(M,X)»RADC4(M,X)+RADZ(M.K) 

182 CONTINUE 
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C TOTAL UPWELLING RADIANCE IN THE SPECTRAL REGION 5-50M 

C 

DO 185 K-16,45 

RADCSCM^-RADCSCM.X^RADZCM.K) 

183 CONTINUE 

RADC5(M,X)-RADC5(M,X)+RADC3(M,X) 

C 

C TOTAL UPWELLING RADIANCE IN THE SPECTRAL REGION 5-200M 

C 

DO 184 K-1 ,195 

RADC6(M,X)»RADC6(M,X)+RADZ(M,K) 

184 CONTINUE 
190 CONTINUE 
C 

C UPWELLING FLUXES AND RADIANCES ARE COMPUTED FOR PARTLY 

C CLOUDY CONDITIONS BY EVALUATING WEIGHED SUMS OF THE 

C CLEAR AND OVERCAST VALUES 

C 

DO 414 NF-1 ,5 

FLUX(NF,X)-(V.-CC(NF))*FLG+CC V 'NF)*FLC 
DO 415 K-1 ,7 

RADCC(M,NF,X)*(1 . - CC ( NF ) ) *RADG (M) +CC ( N? ) *R ADC6 ( M , X ) 

C 

C ANISOTROPIC FUNCTIONS ARE CAUCULATED 

415 AISFUN(M,NF.X)*PI*RABCC(M,NF,X)/FLUX(NF,X) 

414 CONTINUE 
159 CONTINUE 

RDG1 = RDG2-RDG3*RDG4*RDG5*RDG6=»RDG7-0. 

FXG 1 »FXG2«FXG3*FXG4*FXG5=FXG6=FXG7»0. 


C 

C 

DO 170 K-96,195 
rdgi=rdgi+rdg(k) 

FXG1 *FXG 1 +FXG ( K ) 

170 CONTINUE 
C 

C 

DO 171 K*46 , 95 

RDG2=RDG2+RDG(K) 

FXG2*FXG2+FXG(k) 

171 CONTINUE 
RDG3-RDG1 +RDG2 
FXG3 = FXG1 + FXG2 

C 

C 

DO 172 K=76,90 

RDG4-RDG4+RDG(K) 

FXG4=FXG4+FXG(K) 

172 CONTINUE 

C 

c 
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1 73 
C 


174 

C 

C 

67 

68 

69 

70 
76 
64 


65 

71 

7? 


DO 173 K-16,45 
RDG5-RDG5+EDG(K) 
FXG5-FXG5+FXG(K) 
CONTINUE 
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RDG5-RDG5+RDG3 
FXG5-FXG5 + F7.C? 

DO 174 K- 1 ,15 

RDG6*RDG6+RDG(K) 

FXG6»FXG6+FXG(K) 

CONTINUE 

RDG7-RDG6+RDG5 

FXG7-FXG6+FXG5 


WRITE ( 6 , 67) 

FORMAT ( 1 HI /I 9X2 1 HMODEL ATMOSPHERE USED) 

WRITE( 6 , 68) 

FORMAT (///7X3H ALT , 5X5HPRESS , 6X4HTEMP , 6X9HWATER VAP, 8X5 HO ZONE) 
WR:je(6, 69) 

FORMAT (6X4H( KM) ,5X5H(ATM) ,5X5H(KEL) ,9X6H(PPMV) , 7X6H ( PPMV ) // ) 
WRITE(6 , 70) (ALT (L) ,PR(L) ,TE(L) , VMR ( L , 4 ) , VMR ( L , 5 ) ,L-1 ,30) 
FORMAT (FI 0.1 ,F10.5,F10.2,E15.4,E15-4) 

WRI TE( 6 , 76) 

FORMAT (1H1 /25X29HCLEAR ATMOSPHERE PATH LENGTHS///) 

WRITE(J, 64) 

FORMAT (6X2HTH,7X3HULG,8X3HPRG,7X3HTEG, 8X2 HUG , 

1 8X4HPREG, 6X4HTEMG,6X3HALG) 

WRITE (6 f 65) ( (TH(LA) ,ULG(LA,N) ,PRG(LA,N) , TEG(LA,N) , UG(LA,N) , 

1 preg(la,n) ,temg(la,n) , alg ( la , n ) ,LA-1 ,NL) ,N-1 ,5) 

FORMAT (//5(1 6 ( 18 f El 2. 3, FI 0. 5, FI 0. 2 , El 2. 3, FI 0.5, FI 0. 2, F3. 3/)) ) 
WRI TE ( 6 , 71) 

FORMAT (/40X30HINTEGR AT i'J RADIANCE AND FLUXES) 

WRITE (6 , 72) 

FORMAT (/30X "VARIATION WITH NADIR ANGLE (THETA) AND ", "CLOUD 
1 COVER" ) 


WRITE(6 ,84) 

84 F0RMAT(/4OX, "SPECTRAL RANGE -5 - 10 MU") 

WRI TE ( 6 , 73) 

73 FORMAT (///I 8X , 3HRADI ANCE , 4X8HRADI ANCE , 4X8HRADIANCE , 

1 4X8HRADI ANCE, 4X8HR ADI ANCE, 4XSHRADI ANCE, 4X8HRADI ANCE) 

WRITE( 6 , 74) 

74 FORMAT ( 5 X, "CLOUD COVER” , 2X , 8HTHETA+00 , 4X8HTHETA= 1 5 , 4X8HTHET A-30 , 
1 4X8HTHETA-45 .4X8HTHETA-60, 4X8HTHETA-70, 4X8HTHETA-80//) 

WRI TE ( 6 , 60)(CC(X) , (RADC1 ( M , X ) , M- 1 , 7 ) , FLUX (X ) , X- 1 , 2 ) 

WRITE(6 ,71 ) 

WRITE(6 ,72) 

WRITE(6 ,85) 

85 F0RMAT(/4CX, "SPECTRAL RANGE - 10 - 20 MU") 

WRITE(6 ,73) 

WRITE(6,74) 

WRITE (6 ,60) (CC(X) , (RADC2(M,X) , M = 1 , 7 ) , F LUX ( X ) , :l = 1 , 2 ) 
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WRITE(6 ,71 ) 

WRITE(6 ,72) 

WRITE (6,86) 

86 F0RMAT(/4OX, "SPECTRAL RANGE =5-20 MU") 

WRITE (6,73) 

WRITE(6 ,74) 

WRITE (6 ,60) (CC (X) , (RADC3(M,X) ,M=1 ,7) ,FLUX(X) ,X = 1 ,2) 


WRITE (6,71) - 

WRITE(6 ,72) 

WRITE(6 ,87) 

87 F0RMAT(/40X, "SPECTRAL RANGE =5-70 MU") 

WRI1'E(6 ,73) 

WRITE (6 , 74 ) 

WRITE (6 ,60) (CC(X) , (RADC4(M,X) ,M= 1 , 7 ) , FLUX (X ) , X = 1 , 2 ) 


WRITE(6 ,71 ) 

WRITE(6 ,72) 

WRITE(6 ,88) 

88 F0RMAT(/40X, "SPECTRAL RANGE = 10.5 - 12.5 MU") 
WRITE(6 ,73) 

WRITE ( 6 ,74) 

WRITE (6 ,60) (CC(X) , (RADC5(M,X) , M= 1 , 7 ) , FLUX (X ) , X = 1 , 2 ) 


WRITE (6 ,209) EMG , TEMPG , LT , LCB , CEMI ( 1 ) 


WRITE(6 ,71 ) 

WRITE(6 ,72) 

WRI TE ( 6 ,89) 

39 F0RMAT(/40X , "SPECTRAL RANGE = 5 - 200 MU") 

WRITE(6 ,73) 

WRITE ( 6 ,74) 

WRITE (6 ,60) (CC(NF) , ( R ADCC ( M , NF , 1 ) ,M-1 ,7) ,NF»1 , 5 ) 
60 F0RMAT(F12.2,2X,7E12.5/) 

WRITE(6 ,90) 

90 F0RMAT(/44X , "ANISOTROPIC FUNCTIONS"/) 

WRITE (6 ,74) 

WRITE(6 ,60) (CC(NF) , ( AISFUN (M , NF , 1 ) ,M= 1 , 7 ) , NF= 1 , 5 ) 

I F ( LC B . LE.6) GO TO 8585 

WRITE (6 ,209) EMG, TEMPG, LT , LCB, CEMI (2) 

WRITE(6 ,71 ) 

WRITE ( 6 ,72) 

WRITE(6 ,89) 

WRITE ( 6 ,73) 

WRITE(6 ,74) 
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WRITE (6 ,60) (CC(NF) , ( RADCC ( M , NF , 2 ) ,M«1 , 7 ) , NF- 1 , 5 ) 

WRITE (6 ,90) 

WRITE(6,74) 

WRITE (6 ,60)(CC(NF) , ( AISFUN (M , NF , 2 ) ,M-1 , 7 ) .NF- 1 , 5 ) 

8585 WRITE (6 . 78) 

78 FORMAT (‘t HI / I 7X64HCLEAR AND CLOUDY RADIANCES AND FLUXES IN THE 
1 SPECTRAL RANGES///) 

WRI TE (6 , 79) 

79 FORMAT (28X7HRANGE 1 ,8X7HRANGE 2 , 8X7HRANGE 5. 

1 8X7HRANGE 5,8X 1 5HRANGE 4 ( WINDOW) , 5X7HRANGE 6) 

WR ITE ( 6 , 80) 

80 FORMAT (28X7H5_1 0 MU,7X8H10_20 MU,8X7H5_20 MU,8X7H5_50 MU, 
18X14H 10.5- 12.5 MU,5X8H5-200 MU//) 

WRITE (6 ,81 ) RDG1 , RDG2 , RDG3 , RDG5 , RIG4 , RDG7 

81 F0RMAT(7X1 3HCLEAR ATM RAD , 4E1 5 • 5 , E20 . 5 , El 5 • 5 ) 

82 FORMAT (8X1 2H0VERCAST ftAD , 4E 1 5 • 5 , E20 . 5 / ) 

WRITE (6 ,83)FXG1 , FXG2 , FXG3 , FXG5 , FXG4 , FXG7 

83 FORMAT (6X14HCLEAR ATM FLUX , 4E 1 5 . 5 , E20 . 5 , El 5 • 5 / ) 

WRITE (6 , 61) 

61 FORMAT (1 HI /25X45HSPECTRAL TRANSMITT ANCES , RADIANCES AND FLUXE 
WRI TE ( 6 , 63) 

63 FORMAT (///6X4HFREQ,7X8H TRANS T0,6X9HCLEAR ATM , 6X9HCLE AR ATM, 

1 7X8HTRANS TO . 6X9HCLOUD TOP , 6X9HCL0UD TOP) 

WRITE (6 , 75) 

75 FORMAT (6X4HCM_1 , 8X7FSUR FACE, 7X8H RADIANCE, 5X10HFLUX(W/M2), 

1 6X9HCL0UD TCP.7X8HRADIANCE,5X1 OHFLUX ( W/M2 ) // ) 

WRITE (6, 62) (FRC(X) .TRG(K,1 ) ,RDG(K) , FXG(X) ,TRC(K,LC) ,RDC(K) , 

1 FXC(K) , K* 1 , KRT ) 

62 FORMAT (Fi 0.0 , FI 5-5,2E15.5,F15.5,2E15-5/) 

211 CONTINUE 

212 CONTINUE 

213 CONTINUE 
STOP 
END 


SUBROUTINE TRANS(ALB,PL,KR) 

INTEGER W 

COMMON/TRANE/ AVSI (5 , 195), NSI (5 ,195) ,FRC(l 95 ) , TRA( 1 95 ) , 
1X1 (26) ,T1 (26) ,X2(21 ) ,T2(21 ) , DEL A , JD , PI 
C 

C CALUCULATES TRANSMITTANCE IN EACH INTERVAL AT ALL 

C ALTITUDES CONSIDERING DIRECT AND WING CONTRIBUTIONS 

C 
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RHO=ALB/DELA 
DO 300 K-1 , KR 
TRA(K)-1 . 

DO 30 i J-1 , KR 
TRD=1 . 

JA«IABS( J-K) 

if(ja.gt.jd) go to 301 
zi*frc(k)-frc(j) 

EPSI*ZI/DELA 
DO 302 1*1 ,5 
NSJ*NSl(l,J) 

PNUM-RH0*RH0*AVSI(I , J)*PL/(PI*ALB) 
RES*0 . 

IF (J.NE.K) GO TO 303 
DO 304 ¥=1 ,26 

YY«PNUM/(X1 (W)*X1 (W)tRHO*RHO) 

IF (YY.GT.675.) GO TO 305 
Y-EXP(-YY) 

GO TO 304 

305 Y*0 . 

304 RES=RES+Y*T1 (¥) 

GO TO 302 

303 DO 306 ¥=1,21 

YY=PNUM/( (EPSI-X2(¥) ) * ( EPSI-X2 ( ¥ ) ) ) 
IF (YY.GT.675. ) GO TO 307 
Y= EXP (- Y Y) 

GO TO 306 
307 Y*0 . 

306 RES=RES+Y*T2(¥) 

RES*RES/6 . 

302 TRD=TRD*RES**NS J 

301 tra(k)=tra(k)*trd 

300 CONTINUE 
RETURN 
END 


